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PREFACE 


This issue contains the papers* which were presented at the Trans- 
portation Research Symposium held at Purdue University on February 
27-28, 1957. Planned as the first of a series of such symposia (to be 
held from time to time at Purdue University, when warranted by pro- 
gress in transportation research) these papers report problems and 
results from a variety of sources at which relevant research is being 
conducted. These papers, therefore, serve two purposes: they report 
accomplishments to date, along with the status of work being undertaken 
by various groups. They also provide a benchmark for gauging future 
progress in transportation research. 


The papers in this issue of the Naval Research Logistics Quarterly 
are presented in the same order used for the Symposium. Commencing 
with an historical summary by Roger Crane, the series concludes with 
a report of a new method of computation for general distribution prob- 
lems by William Prager. The logic and balance of the Symposium’s 
organization is thereby preserved, except forthe record of participants’ 
remarks and the concluding statement by Professor A. Charnes, Chair- 
man of the Symposium. Because of their wealth and variety, it was not 
possible to condense these remarks sufficiently for presentation in this 
issue of the Naval Research Logistics Quarterly. 








Acknowledgment is due the Institute of Management Sciences, the 


Operations Research Society of America, and the American Institute of 
Industrial Engineers, Inc. These professional societies served as 
sponsors of this Symposium conducted by the Management Sciences 
Research Group of the Department of Industrial Management and Trans- 
portation and the Division of Adult Education of Purdue University. 


Acknowledgment for assistance in preparing these materials for 
publication is due the Management Sciences Research Group at Purdue 
University, and especially to the director of this group, Professor 
Abraham Charnes. 


M.I. R. 





*Since their report was not available, the paper “Communication — 
Transportation Networks” by Kalaba and Juncosa has been omitted 
from this series in order to avoid further delay in publication. 








SOME RECENT DEVELOPMENTS IN TRANSPORTATION RESEARCH 


Roger R. Crane* 





| A survey of recent history of and techniques in transportation methods. | 





INTRODUCTION 

It is certainly a pleasure to be here and to have the opportunity to address the first 
symposium on transportation research held at Purdue University. Purdue is justly renowned 
for its pioneering work in transportation research and engineering, and it is entirely fitting 
that these symposia be commenced at this university. 

This is the first meeting of this symposium, and because our audience consists of both 
scientists and businessmen, I have selected for this presentation the topic: Recent Develop- 
ments in Transportation Research. In this presentation we will start by discussing the nature 
and importance of transportation research and follow this with a brief description of important 
recent developments in this field. 


THE NATURE AND IMPORTANCE OF TRANSPORTATION RESEARCH 

When someone mentions scientific research to most business people, it has been my 
experience that they immediately think of the technical research and development leading to 
new and improved products. By product they usually mean some device or form of hardware. 
For many companies this is, indeed, their primary product. It is this product, or these prod- 
ucts, which must be improved and which must be added to if the company is to maintain its 
position in the sales market and to continue to maintain a sound financial position. There are, 
however, many companies in industries, such as transportation, where the primary product is 
not a piece of hardware at all, but a service. Research—scientific research—on this service 
is as important to the financial health and economic survival of these companies as technical 
research and development is to manufacturing. On the face of it, this observation hardly seems 
acontroversial one. Yet it is only in very recent years that this type of "transportation 
research" has begun to be emphasized. 

The meaning of research in this context is clear, but what do we mean by "transporta- 
tion"? We can look around the audience here and see people who have come to this meeting as 
representatives of airlines, pipelines, bus companies, railroads, and the like. Certainly all of 
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you consider that your organizations are in transportation. It is important from the point of 
view of the research man, as we will point out later, that we also understand transportation in 
a somewhat broader context than might be concluded from a review of corporations whichform 
the transportation industry in the United States. We mean by "transportation" the process of 
moving men and material from one location to another for purpose of increasing the utility of 
the object moved, Thus, it is natural to include as "transportation" all forms of materials 
handling and personnel movement, even that carried out, say, within a manufacturing plant. 
The value of emphasizing this generalization is to underscore the similarities among these 
different forms of transportation so that we will keep in mind that research advances in inven- 
tory control, for example, may be applicable directly to the problem of freight car utilization, 


RECENT DEVELOPMENTS 

Certainly transportation research is not new. Scientific methods have been used in 
transportation problems, as in many business problems, for a long time. But recently there 
has been increasing emphasis on this work. One of the most important indications of this is 
the rapid growth of industrial operations research since World War II. During World War II, 
operations research resulted in introducing into the study of business problems many com- 
petent scientists. These scientists were brought into direct contact with operational problems 
during the war. As was only natural, when World War II ended they reasoned that the same 
skills and approach which had been used so successfully in solving military problems might 
well contribute to the solution of the operational problems of American business. 

The growth of these new techniques and this new approach in business was slow. It was 
not until about 1950 that any sizeable, organized effort in this area commenced. The growth of 
OR in industry since 1950, however, has been limited primarily by lack of qualified personnel. 
Since 1950, two societies, devoted primarily to the application of scientific methods to business 
problems, have been founded and have successfully grown to considerable size. These are the 
Institute of Management Sciences, which was founded in 1953, and the Operations Research 
Society of America, founded in 1951. Both have memberships well in excess of 1000, and The 
Institute of Management Sciences has many additional members abroad. Both publish Journals, 
one called Management Science and the other Operations Research. You may be interested to 
know that two of the people most influential in the establishment of The Institute of Management 
Sciences are sitting here at this table today - Professor Charnes, of Purdue, on my left, and 
Professor Cooper, of Carnegie Institute of Technology, on my right. 

When I was asked several months ago by Professor Charnes to undertake the task of 
addressing you at this opening luncheon of the first Transportation Symposium, I felt reason- 
ably well equipped to do so, having been actively engaged in Operations Research on trans- 
portation problems of several years before joining the Public Accounting firm of Touche, 
Niven, Bailey & Smart at the end of 1955. Until 3 months ago, I had been somewhat preoc- 
cupied, as you might well imagine, becoming familiar with public accounting and had little 
time to keep abreast of new developments in transportation research. I am pleased to report 
that during the past three months I have again had the opportunity of working directly on 
research projects for two transportation systems—an airline and a railroad—and am also 
connected with a project for a bus company. 

This research work, plus the study brought about by the request of Dr. Charnes, have 
led me to a surprising realization. As I moved back into the field to see what had developed 
since 1955, the ‘mpact of the increase in the rate of scientific research on transportation 
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problems seemed to be tremendous. It was hard to believe that I had been away from this 

field for only about one year. I reviewed the literature which was available, primarily the 
journals of the Institute of Management Sciences and Operations Research Society of America, 
and talked to many of my friends in transportation research about new developments in their 
organizations. To obtain as complete and authoritative a picture as possible of these develop- 
ments and others since 1950, I wrote to about a dozen organizations which, to my knowledge, 
were engaged in Operations Research as applied to transportation. Each of these individuals 
were asked for several paragraphs summarizing their efforts and listing some of the types of 
problems in which they had been engaged. Their replies form a basis contributing significantly 
to the following: 

In 1951 the Case Institute of Technology established its Operations Research group as 
part of its Engineering Administration Department. Though it is not generally known, this 
group was started largely through the efforts of Mr. John Kusik, financial vice president of 
the Chesapeake & Ohio Railroad. Mr. Kusik interested Case Institute in acquiring certain 
qualified individuals desirous of applying the scientific method to business problems, As a 
result of Mr. Kusik's discussions with Case Institute, Professors C. West Churchman and 
Russell Ackoff joined the Case staff and commenced their now well-known Operations Research 
Group. 

Another point which may not be too well known, is that the first problem on which Case 
Institute worked was a joint project with the Chesapeake & Ohio Railroad. The subject of this 
study was the application of statistical sampling to interline settlement of freight revenue. 

Since 1951, when it successfully commenced its operations research work with the 
Chesapeake & Ohio Railroad, Case Institute's Operations Research Group has done research 
for a number of transportation concerns, including the Railway Express Agency, a major 
domestic Airline (to be reported by Dr. McCloskey of Case Institute this afternoon), and an 
analysis of the distribution of raw materials for a major metals processor. As pointed out 
earlier, such studies as the latter, although made for a manufacturing company, can very often 
contribute significantly to improved transportation service itself. 

The financial department of the Chesapeake & Ohio Railway, which had commenced its 
research with Case Institute, continued in this field, adding a number of highly skilled person- 
nel to its staff and widening the scope of its research to include studies of information and data 
handling for the entire system and the analysis of potential applications and use of large-scale 
electronic computers. All of you are probably aware that the C. & O. has already installed a 
large-scale computer which they are using not only for financial operations but also for certain 
operating problems of the railroad. 

Almost entirely apart from the development of Operations Research in the financial 
department of the C. & O. just described, the maintenance of way and structure department of 
the C. & O. joined together with Melpar, Inc., to carry out a number of operations research 
studies on rail failures and maintenance. They went beyond this to detailed studies of other 
maintenance problems, particularly the optimum utilization of work equipment and manpower. 
The C, & O. now has a number of members of their staff who are experienced in and capable 
of carrying out operations research work on maintenance problems. 

In 1951 the University of California at Los Angeles established its Institute of Trans- 
portation and Traffic Engineering. The institute has since carried out considerable research 
o highway and sea traffic problems by the aid of computers and on the development of mathe- 
matical models of highway and street traffic flow. In 1951 UCLA also started its cargo handling 
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research project sponsored by the Office of Naval Research and the Maritime Administration, 
The purpose of this continuing research project was the study of all forms of transportation 
of material by land and by sea. Mr. Russell O'Neill, who was the first project director of this 
group, continues in this capacity and will appear on this program at Purdue. 

In 1952 Melpar, Inc., the research subsidiary of Westinghouse Air Brake Company, 
started an operations research group with the objective of studying all forms of transportation 
and, in particular, the railroads of the U.S. I was connected with this group almost from its 
inception. By the middle of 1953 this group had grown to 12 members and was carrying out 
operations research projects with five major U.S. railroads in such diverse fields as maintenance of 
way and rolling stock, classification yards, signalling, and terminal area operations. Some 
studies also were made of other forms of transportation, such as highways and streets and 
aircraft. 

In 1954 Westinghouse Air Brake Company decided to decentralize its research activi- 
ties, including its operations research work. To carry on a longer range aspect of the projects 
developed by the Melpar Operations Research Group, a grant of $75,000 was given by Westing- 
house Air Brake Co, to the Graduate School of Industrial Administration of Carnegie Institute 
of Technology. You might almost say that if this grant had not been made, you would not be 
here. I say this because I am sure that it is, at least in part, the result of this grant that 
Professor Charnes, now at Purdue University but then at Carnegie Tech, became interested 
in transportation research. Professor Charnes made a number of contributions to the research 
work carried out on transportation at Carnegie Tech. 

Carnegie Tech established a Transportation Research Group to study transportation 
problems in general and, like the Melpar Research Group, to look in particular at railroad 
problems. This group is actively engaged in a number of research projects. Professor Harold 
Wein has told me something of the projects which are under way at Carnegie Tech and promises 
that these results will be published soon. The work of his group includes studies on the sched- 
uling of engine crews, costs in classification yards, yard location, and the measurement of 
productivity and efficiency in railroad transportation. Professor Miller of Carnegie Tech and 
Mr. Soyster of the Union Railroad in Pittsburgh will be talking tomorrow about a part of this 
research. 

One of the most important recent developments was the great interest which, since 
1953, the young Railway Systems and Procedures Association began to show in operations 
research. At their meetings during and since 1953, the RSPA has sponsored many presenta- 
tions by operations research scientists. The RSPA started out primarily in the financial and 
accounting areas but has expanded its interests and participation to include all important 
aspects of railroad operations, including traffic, maintenance, and long-range planning. Heavy 
emphasis has been placed upon the role of the electronic computer in railroad operation. A 
former president of the RSPA, who active during this early period, is here today. He is 
Mr. B. E. Wynne, controller of the Western Maryland Railroad, 

In 1954 transportation research activity continued to increase. United Airlines had long 
been interested in the application of scientific methods to their operating problems. Under the 
leadership of Warren Alberts, who will be chairman of this afternoon's session, it began its 
first full-scale operations research project. The purpose of this project was to study aircraft 
routing and line maintenance; but as the project continued, the research group began to look at 
virtually all aspects of the airline's operations. At this point they have developed a model 
which dynamically simulates all maintenance and operations at various stations. Quality and 
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service levels as related to manpower control have been investigated, as have inventory control, 
reservations, forecasting, manpower control, and accounting. 

Midwest Research Institute and Stanford Research Institute, both intensified their trans- 
portation research in 1954. Midwest Research Institute worked with the Great Lakes Pipeline 
Co, on scheduling the flow of products through a pipeline and analyzing congestion in terminal 
areas. Midwest Research Institute has studied the utilization of tractors and semitrailers for 
atrucking company and the reservation office service of TransWorld Airways. They have also 
done research on distribution problems for International Harvester. 

In 1954 Stanford Research Institute, whose industrial economics group long had been 
involved in studies of transportation problems, began a major operations research project for 
the Southern Pacific Railroad. This project was at first primarily in inventory control for the 
purchases and stores department. Several reports have been written on the results of this 
project and are available in the OR literature [29]. More recently Stanford Research Institute 
and Southern Pacific have been focussing their attention on the fundamental transportation 
problem—the utilization of equipment, engines, and rolling stock. 

This latter problem—that of getting the most return on the heavy investment in vehicles, 
motive power, and crews—is common to most forms of transportation. It has bothered the 
railroads in particular, because their freight cars move loaded less than 10 percent of the 
time. Mr. Feeney, of Stanford Research Institute, has developed a mathematical approach to 
this problem and a procedure for applying this approach, by use of a medium-sized computer, 
to improve the utilization of empty freight cars. Mr. Feeney reported this new development 
only last month at a Case Institute of Technology seminar. 

Stanford Research Institute has analyzed a number of distribution and logistics prob- 
lems for the military and, in the process, has developed ideas which ultimately may be of 
considerable value to American industry. 

In 1954 Leslie Edie, head of transportation research at the New York Port Authority, 
published his now famous paper on traffic delays at toll booths [6]. This paper was awarded 
the first Lanchester prize of the Operations Research Society of America, a yearly prize 
awarded to the outstanding operations research paper published. The screening is an exhaus- 
tive one, and the fact that this award was made to Mr. Edie for a transportation application is 
a significant indication of the interest in, and the potential of, transportation research. In his 
work Mr, Edie has made particularly good use of queueing or waiting-line theory, a mathemati- 
cal theory of congestion of items such as automobiles waiting to be serviced at the toll booths, 
which relates the congestion to the facilities available for the servicing of the items. 

Still in 1954, the Bessemer and Lake Erie Railroad, which, as already mentioned, is 
closely connected with the Union Railroad of Pittsburgh, became involved in two research 
projects: the determination of the location of major classification yards and the analysis of 
main-line shipments. 

In 1955, under the leadership of its financial vice president, Mr. Phillip Small, Pacific 
Inter-Mountain Express entered into the operations research field; first, by entering into a 
project with a consulting organization and, second, by hiring an operations research scientist 
from Case Institute. This research man has since been involved in a variety of operational, 

48 well as financial, problems. 

Westinghouse Air Brake Company carried out an operations research project on the 

Scheduling of barges and tugs for a major barge line. 
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The Western Maryland Railroad acquired Mr. Wynne from the Bessemer and Lake Erie 
Railroad in Pittsburgh. His familiarity with operations research soon led to the initiation of a 
full-scale OR project at the Western Maryland. The major objective of this project was to 
obtain maximum utilization of the large and expensive classification yards scattered about his 
system. 

In 1955 Professor Charnes came to Purdue University and started his Management 
Sciences Research group with which you are all familiar. This group is active in transporta- 
tion research. Professor Bartlett of this group and Mr. Gleaves of the Frisco Railroad will 
talk this afternoon about one of its major transportation projects. 

Perhaps the most important development in 1956 was not the completion of a successful 
research project at all, but the recognition of the value of scientific research on transportation 
problems to a transportation company. This has been demonstrated by the appointment of two 
well-known operations research scientists to outstanding administrative posts in transportation 
systems. Dr. Omand Solandt, formerly in military operations research in England and Canada, 
was appointed assistant vice president of research and development for the Canadian National 
Railways, and Dr. Foster Weldon, of the Operations Research Office of the army, who will be 
chairman of tomorrow morning's session, was appointed director of research for the Matson 
Steamship Lines in San Francisco. Both of these men were selected, in part, because of their 
demonstrated ability to apply the scientific method to operational problems. I believe this 
indicates a significant new trend in the transportation industry, perhaps the growth of trans- 
portation research departments of scope and stature comparable with our great industrial 
technical research and development departments, 

Unfortunately, in this brief presentation I have had to leave out many organizations that 
are active in the transportation field. Although the list could never be complete, a few addi- 
tional organizations should be mentioned in passing. First of all, we should mention the Rand 
Corporation in Santa Monica, where several hundred scientists are engaged, in part, in the 
study of various forms of transportation, including railroad, shipping, and aircraft. 

We should mention Brown University, where Dr. Prager, who will address us tomorrow, is 
located, and the Arthur D. Little Co., which has carried out a number of transportation research 
projects. 


PUBLICATIONS 

One of the earliest papers published on transportation in the Journal of the Operations 
Research Society was John Everett's article on the application of queueing theory of overseas 
shipping. Although queueing theory or waiting-line theory had long been known by scientists, 
Everett's application to shipping opened a new field for the operations research scientist. 
Notable among subsequent papers were those by Charnes and Cooper and Dantzig, which develop 
linear programming and its application to distribution and transportation problems. One of the 
special techniques connected with linear programming is, in fact, called "the transportation 
method," 

I mentioned the RSPA earlier. There are a number of articles on operations research 
in the proceedings of their meetings. The Railway Systems Procedures Association had a 
seminar with Railroad officials, which was specifically devoted to OR in that field. The 
American Railway Engineering Association has been studying OR for some time and recently 
published in their Bulletin an article [27] on this subject as it pertains to transportation. Both 
TIMS and ORSA have devoted sessions at their annual meetings to transportation. Publications 
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on transportation will undoubtedly continue to increase rapidly as interest in the field increases, 
Abrief bibliography is appended, 


SUMMARY 

Perhaps we can summarize this presentation by listing the types of transportation 
industries which have been involved in the use of these scientific methods in the solution of 
their business problems, the problem areas in which these studies have been made, and the 
techniques which have been used. To some of you, a mere listing of techniques may be of 
little value, though it may help to give you all a better idea of the interrelationship between the 
techniques and the problem areas and of the tremendous potential for this type of work in trans- 
portation. 

The following activities already have been involved in scientific research on transpor- 
tation problems: airlines, railroads, barge lines, pipe lines, truck lines, and bus lines—vir- 
tually all kinds of transportation. The following problem areas have been attacked: Congestion, 
scheduling, reservations, paper work, maintenance-of-way and equipment, accounting, facilities 
planning, and personnel. The following techniques have been used: linear programming in 
several forms, trans-shipment and the assignment technique, the travelling-salesman technique, 
nonlinear and dynamic programming, waiting-line theory, and the Monte Carlo method. In fact, 
most major branches of mathematics, as well as electronic data-processing, have contributed 
something to transportation research. 


CONCLUSION 

We would like to conclude with these three points: First, that transportation is the 
physical movement, rapidly and safely, of men and material from one place to another. These 
movements can be measured and hence form an excellent basis for the development of scien- 
tific theories. Second, service is what transportation systems sell. Research on this service 
is research on the basic product sold and can lead to new and improved products and hence to 
improved financial health of a transportation company. Third, all forms of transportation, 
from a research point of view, have much in common; successful research in one area can 
frequently contribute materially to successful research in another superficially unrelated 
area, For this reason progress in transportation research may well be rapid; but only if 
symposia, such as this one, are held to bridge the communications gap between the diverse 
types of transportation represented here. Thank you. 
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AN ANALYSIS OF STEWARDESS REQUIREMENTS AND SCHEDULING 
FOR A MAJOR DOMESTIC AIRLINE 


Joseph F. McCloskey and Fred Hanssmann 


Case Institute of Technology* 





This paper studies the optimal utilization of stewardesses and optimal 
allocation of trips to bases by a special approximation to the solution of 
a linear programming problem. 











This is a summary of an operations research study carried out by the Operations 
Research Group of Case Institute of Technology for one of the leading domestic airlines of the 
United States. As is the case with practically all of the OR work undertaken by Case for indus- 
try, this study had the mixed objectives of acquainting the airline with possible uses of OR; of 
training selected personnel of the company in the methods, techniques, and tools of operations 


research; and of providing the company with solutions to one or more problems. 

The work covered a period of about a year, kept an average of three Case and two air- 
line staff members occupied, and was concerned with three interrelated aspects of the airline's 
operations. In general, the three objectives of the study were met. From the OR point of view, 
the two most interesting aspects of the study were the technique developed for handling the 
assignment of trips to bases and the demonstration of the value of having flexibility within an 
OR group. 

As will be seen, the three parts of the problem each required a different set of tech- 
niques for solution. The first involved some rather complex differential calculus; the second, 
a considerable variety of statistical techniques; and the third, linear programming. This 
necessitated more changes among the junior members of the team than is usually desirable, 
but it does support the argument of those who hold that OR should be done by specialists with 
deep knowledge of a particular field. Fortunately, we had the specialists, but the excessive 
turnover militated against the development of an overall understanding of the system by all of 
the team members. This, in turn, made the work go more slowly than either we or the airline 
liked. The airline people responsible for the study are to be congratulated for their patience 
(not always evident when delays were encountered); thanks to that patience, all concerned were 
Smiling and happy after the final briefing on the final report. 





*In a genuine sense, this paper is a group effort. Special mention must be made of the work 
done by Eliezer Naddor, now of The Johns Hopkins University 
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We also wish to thank the airline for permission to reveal certain information neces- 
sary to an understanding of the techniques used in the study. All numbers have been altered 
in order to protect the client. 


THE SCHOOL-SCHEDULING PROBLEM 

The initial problem presented to the research team by the airline was that of scheduling 
the stewardess school operated by the airline. Here, the capacity of the school; the relatively 
long lead time involved in selecting, processing, and training the stewardesses; the seasonal 
fluctuation in requirements for stewardesses; and the high turnover rate were the major com- 
plicating factors, Illness, vacations, leaves of absence, the preference for June and September 
(both bad in terms of seasonal fluctuations) as months in which to get married, and the need 
for special training (as in making the transition from one aircraft type to another) all had to be 
taken into account, as well. 

The analogy between this type of personnel problem and the "inventory" problem 
quickly became apparent, but this presented the additional problem of determining the costs 
of overages and shortages. This, in turn, brought out some unusual aspects of the system that 
might otherwise have been overlooked but which led rather directly into other phases of this 
study. 

Thus, it was difficult to assess the cost of a shortage, because no flight had ever been 
cancelled for lack of a stewardess. Numerous courses of action, including temporary re- 
employment of a "retired'' stewardess, are available to a crew scheduler when a shortage 
threatens; it was through an examination of the costs of some of these courses of action that a 


suitable cost of shortage was determined. On the surplus side, the system can absorb rather 
large numbers—a few at each base—in useful work. When the surplus problem becomes acute, 
careful (and acceptable) manipulation of vacations and leaves of absence can take care of a 
large part of the surplus. This flexibility of the system made it possible to limit consideration 
of costs to a fairly narrow range of possibilities. Within this range, it was found that the cost 
of surplus was very close to wages paid and that the cost of shortage was roughly three halves 
of this amount. The problem may then be conceptualized as in Figure 1. 














—— Lead Time Period»}!~—__———_ Planning Period ——> 


Planning Graduation 
Time Date 


S = level at planning date X = losses in lead time period 
N = additions in lead time period Y = losses in planning period 
n = number in class in question z = auxiliary variables S -R + N + n) 


Figure 1 - Stewardess level during planning period 
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The mathematical model and its solution are discussed in the first annex to this paper. 
Inasmuch as there was little probability of finding the mathematical capability involved in 
handling this model among the personnel responsible for the school, it was necessary to devise 
a practical operating technique. The solution was to devise a table of z values (see Fig. 1) 
such that the operator could obtain the appropriate z value by reading across to the expected 
loss in the planning period and then down to the expected loss in the lead period. The z value 
thus obtained could then be used to solve the equation, 


n=z+R-S-N, 


where n is the optimal class size, R is the requirement for stewardesses, S is the number 
available, and N is the expected graduations. 

An example will make clear the use of this equation. Suppose that planning is being 
done on December 31 for a class expected to graduate from the school on February 28. Gradu- 
ation will produce 35 stewardesses on January 2 and another 23 on February 1. There are 965 
available stewardesses and the requirement for March has been set at 975. Expected losses 
(which are known within rather narrow limits for the time intervals involved) are 38, 42, and 
48, respectively, for the first three months of the new year. Thus, expected losses are 80 in 
the lead period and 48 in the planning period. The interpolated value from the table for these 
two figures is 110. Therefore, 


n=Z+R-S-N 
n= 110 + 975 - 985 - 58 
n= 42. 


Under some circumstances, n may be negative, or too small to justify a class; in such 
cases, management has the option of foregoing one or more classes until the value of n justifies 
resumption. Under other circumstances, n may be so large as to exceed the capacity of the 
school; in such cases, the school must be operated at capacity until the value of n drops below 
capacity. (In practice, application of this equation has been extended to anticipate major sea- 
sonal fluctuations so that very small and very large n's are avoided. A simple decision rule 
also has been developed to determine when one or two classes should be scheduled and to 
anticipate when a third class may be required in the school at any one time.) 

With a procedure thus available for determining the desired size for any one class, the 
next problem is that of class frequency. Analysis showed that the school is so organized that 
the incremental costs of having two classes in the school at one time are negligible unless the 
Capacity of the dormitory is exceeded by a considerable amount; additional facilities would be 
required to handle a third class of reasonable size. This led to the conclusion that classes 
should be scheduled as frequently as two classes could be accommodated. In the interest of 
controlling shortages and surpluses, the classes should be spaced as evenly as possible through 
the year, but practical considerations, like shutting down the school over the Christmas holi- 


days, require that the second class enter the school before the first has reached the halfway 
mark, 


In order to get some measure of the effect of class frequency on the system, tests were 
Tun on the school output over several previous years. These tests showed that 20 classes per 
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year reduced the costs of surplus and shortage 67 percent below those for 10 classes per year, 
Thus, given the present organization of the school and the same output of graduates, more 
frequent classes provide finer control of surpluses and shortages and, because of smaller 
class sizes, make possible more personalized and effective instruction. 


DETERMINATION OF REQUIREMENTS FOR STEWARDESSES 

An examination of the equation developed for determination of the size of classes in the 
stewardess school shows that all of the inputs but one are known within rather narrow limits, 
The one exception is the system requirement. Determination of the requirements for stew- 
ardesses became, then, the second part of this study. 

As the system operates, flight schedules are prepared, the flights are then assigned to 
bases in the system, and the assigned flights are manned out of each base. Once the number of 
stewardesses required to man the flights out of a base is known, a safety factor—in the form of 
reserve stewardesses—is added. In general, it is expected that the reserves will fly consider- 
ably fewer than the 85 hours per month that they may fly legally. 

There are several alternatives open to the crew scheduler in making the necessary day- 
to-day adjustments in the crew schedules, Thus, one policy might be to reschedule a regular 
stewardess every time she became underprojected (missed a flight because of illness, cancel- 
lation, etc.). Another might be to schedule each regular stewardess for more hours than she 
could fly legally, relying on reserves to take over the surplus in the event there were no can- 
cellations or illness during the month. (In good-weather months this could create serious peak 
requirements for reserves at the end of the month and so would require a much more compli- 
cated system of scheduling.) Another might be to use reserves at every opportunity, up to the 
legal limit of flying hours, before calling on underprojected regular stewardesses at all. 
Another might be to adjust a stewardess' pay to the number of hours flown, thus encouraging 
regular stewardesses to go on the reserve list when they became underprojected. 

"Paper" experiments were devised and run to test several of these alternatives in order 
to get some measure of the reliability of the estimates of system requirements and their sen- 
sitivity to the various alternatives. As might be expected, rescheduling of underprojected 
regular stewardesses led to a reduction in reserve requirements but would have been difficult 
to implement because, as we shall see, stewardesses prize their days off and don't like to have 
their plans upset by being placed in reserve. And, as anticipated, overscheduling created 
special problems. The other alternatives proved less efficient or involved contract negotiations 
that would have deferred implementation. 

Because the airline prides itself on the working conditions that it affords its stewardesses, 
an effort was then made to keep the present conditions intact and to see whether or not the 
reserve requirements could be reduced. This approach had the obvious advantage of requiring 
minimum retraining of base personnel and minimum effect on the day-to-day lives of the 
stewardesses, 

The first approach was to test the validity of the established reserve requirement. This 
was found to fluctuate rather considerably; that is, the ratio of actual to scheduled flying by 
regular stewardesses varied to an extent that the established safety margin appeared neces- 
sary. The sensitivity of this ratio to actual versus planned flights was then tested for all 
months of the year and all bases in the system. It was found that the ratio of actual to planned 
flights varies from slightly over 1.0 to slightly under 0.9, with the better -than-"perfect" per- 
formance coming in good-weather, high-traffic months when extra sectigns of scheduled flights 
are flown. The poorer performance is, of course, almost wholly a function of bad weather. 
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The interesting phenomenon discovered was that when these two computations are 
combined (that is, when the quotient of the actual to the planned flying hours of the steward- 
esses is divided by the quotient of the actual to the planned flights) a very nearly constant 
ratio is obtained. In other words, regardless of weather or base, regular stewardesses 
"capture" a very stable percentage of the actual flying hours, leaving the remainder to the 
reserves, At first glance, this appeared to confirm the validity of the established reserve 
figure. The fact is, however, that only about 55 or 60 percent of the regular stewardesses are 
on duty on any given day (inasmuch as a working month consists of about 13 to 18 duty days), 
so that, in effect, the established reserve was much closer to a reserve about two-thirds again 
as large. 

Thus, by using the stable "capture" percentage in planning, the expected utilization of 
reserves may be increased significantly and the requirement may be reduced accordingly. The 
precise reduction is difficult to determine because of the dynamics of the system and the exist- 
ence of bases from which relatively few flights are flown, but the recommendation for the use 
of a higher planning figure for reserve utilization has the advantage that it may be approached 
incrementally until difficulties are encountered. 


ASSIGNMENT OF FLIGHTS TO BASES 

Simultaneously with the effort to reduce reserve requirements, research was carried 
out on the problems associated with optimum utilization of all stewardesses. Preliminary work 
reduced this to the problem of optimal assignment of flights to bases, The aircraft routings and 
the airline schedule of flights were taken as "given," so that the only variable that could be 


manipulated was time on the ground away from the locations at which the stewardesses were 
based (hereafter referred to as 'away-from-home time"), Minimization of this time has the 
dual advantages of making stewardesses more available for emergency or routine make-up 
flying and of reducing the airline's outlays for meals and lodging of personnel away from their 
home bases, Significant reductions in this away-from-home time can, of course, result in 
some reduction in the number of required stewardesses., 

A more straightforward way of handling this problem would have been to use away- 
from-home costs directly, but this would have entailed considerable delay, inasmuch as crew- 
expense accounts combine in-flight and on-ground allowances. Accordingly, away-from-home 
time was used and was then converted into costs on a sample basis. Use of away-from-home 
time had another advantage from the point of view of implementation: An analysis of the pref- 
erences of the stewardesses themselves showed that the only factor that influenced their bidding 
significantly was the number of free days, at the home base, that the bid afforded. 

A simple illustration will serve to demonstrate the fundamental nature of the problem. 
Consider a flight system that involves only three locations on a line and three flights in each 
direction, as shown in Figure 2. 

If turnaround in either direction at L, is excluded, the possible combinations of the 12 
flights may be represented by a table of the type shown in Figure 3. In this matrix, the cell 
marked ''x"' represents the combination of the four flights 1,2,A,a. If this matrix is to be used 
48a means of minimizing away-from-home time, the away-from-home time associated with 
tach possible combination must be inserted in each of the cells. Then, those three cells must 
be selected that include each of the 12 flight segments (once, and only once), and the sum of 
Whose away-from-home times is minimum. For a small problem of the sort used in this 
illustration, the answer may be obtained by trial and error; but for a realistic problem such as 
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Figure 2 - Simplified flight pattern 
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Figure 3 - Matrix of possible flight combinations 


that represented by the airlines system, trial and error could not yield a minimum solution in 
any reasonable amount of time. And there is no ready-made mathematical technique available 
for solving the problem. Therefore, it was necessary to use an existing mathematical proce- 
dure which yields a good approximation of the minimum and which can be adapted to rapid 
computation of the actual problem. (Rapidity is not always desirable where a research solution 
is involved, but it becomes necessary where the procedure is to be incorporated into a highly 
dynamic operating situation.) 

The method of solution adopted (the linear programming solution of the assignment 
problem, discussed in Annex B) handles pairs of flights or flight segments between two stopping 
points, whether or not these stopping points are bases in the system, Consequently, the prob- 
lem of flight pairings and allocations is reduced to a series of two-city problems, for each of 
which minimum away-from-home time can be determined, If the three non-stop flights in each 
direction between B and A are taken as an example, the logic of the procedure that was even- 
tually developed becomes clear, 

Figure 4 lists the three flights in each direction, with their departure and arrival times. 
The matrix in Figure 5a shows the layover times at A for each pair of arrivals and departures. 
For ease of computation, these times are given in minutes (the upper number in each cell), but 
for ease in relating layover times to restrictions on crews, they are also given in hours and 
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minutes (the lower number). The matrix shows, for example, that a crew arriving on Flight 1 
and departing on Flight 2 will have a layover of 17 hours and 15 minutes. Figure 5b gives the 
same information for B. For the same pairing of Flights 1 and 2, the layover time is 14 hours 
and 45 minutes. Considering only this pairing, the crews involved should be based at B, 
pecause the layover time there must be increased by 24 hours in order to meet legal require- 
ments, while the layover at A is sufficient, as it stands, to meet legal requirements. 
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Figure 4 - Schedule between A and B 
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Figure 5a - Layover time at A 
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The "Third Matrix" 


Figure 5b - Layover time at B 
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When each possible pairing is subjected to this type of step-by-step analysis, the "third 
matrix," shown in Figure 5c, results. There, for the pairing of Flights 1 and 2, the layover 
time at A is inserted and the notation is added that the crews should be based at B. This 
"third matrix" is solved by selecting one cell in each column and each row such that the sum 
of the layover time is minimum. The cells marked with an asterisk (*) in Figure 5c satisfy 
these requirements. This solution is compared with the actual airlines trip allocations in 
Table 1. The computed minimum (the OR Trip Allocation) is seen to be 870 minutes per day 
less than the actual trip allocation. 


TABLE 1 
Layover Times 
(Actual Trip Allocations vs. Computed 
OR Trip Allocations 





Layover Time 





(min) | (hr + min) 





Actual Trip 2 1245 20:45 

Allocations 1755 29:15 
1755 29:15 
Totals 4755 79:15 








OR Trip 5 | 1005 16:45 

Allocations 1485 24:45 
1395 23:15 
Totals 3885 64:45 














Difference 
Saved Daily 870 14:30 

















It is apparent from this example that the mathematical technique may be applied only 
when the matrix is square, because of the requirement that one cell be selected in each column 
and row. Therefore, before the procedure could be developed for the whole system, all flights 
and flight segments had to be "squared" between each pair of cities so that the number of 
arrivals and departures from each city is equal to the arrivals and departures from the paired 
city, even though the routings or number of segments may differ between the two members of 
the pair. Thus, if there are three direct flights from A to B, but four from B to A, it is neces- 
sary to put together other (legal) flights or segments, the first leg of which originates in A and 
the last of which terminates in B, in order to "square" the third matrix for A and B. In other 
situations, it may be necessary to break a current flight into two or more segments in order to 
accomplish the squaring of the matrices. 

Obviously, considerable subjectivity is introduced by this procedure. It is important, 
therefore, that the matrices be squared by personnel thoroughly familiar with the system 
before they are incorporated into the operating system. 
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When this squaring had been accomplished, the procedure was applied to the daily flight 
schedule for November 1955 by using hand computing methods. The flight schedule, the equip- 
ment routings, and the present bases of the airline were held constant in the first trial; apart 
from the squaring, oly the assignment of flights to bases was altered to achieve the desired 
minimum layover time. This hand computation showed that reductions of as much as 20 per- 
cent in away-from-home time for stewardesses might be possible; however, no definite con- 
clusions could be drawn from this trial run, because only the most important restrictions 
could be incorporated in the time available. 

’ A more elaborate model was then developed to include all of the restrictions that 
affected away-from-home time, and a later schedule was then programmed into an IBM 650 
computer. The "squaring" had to be done in advance, but the computer could handle the tasks 
of scanning the first and second matrices, applying the restrictions, and producing the third 
matrix. The third matrices could then be "solved" by the assignment-problem technique or by 
the "map-tack method" as described in the second annex to this paper. Over 100 matrices, 
ranging in size from 1 x 1 to 30 x 30 were involved. 

The additional restrictions reduced the potential saving in away-from-home time 
drastically, but the money involved is such that even the remaining modest saving was worth- 
while. The program was rerun several times, with the restrictions relaxed one at a time, in 
order to provide some idea of the relative costs of the various restrictions. It was also rerun 
with additional bases inserted into the system in order to show what additional savings might 
be possible. This rerun showed that two additional cities served by the airline should be con- 
sidered as bases, given the existing pattern of flights. (As a matter of fact, another base may 
be added to the system—in neither of the recommended cities—because of additional routes 
since granted to the airline.) 

One of the most interesting results of this machine analysis was the discovery that the 
present system of assignment, largely empirical, produced the optimum away-from-home time 
for two-engine flights and that all of the potential saving was in the four-engine assignments. 
No satisfactory reason for this difference was discovered, but there is a possibility that it was 
associated with length of flight or with the method of inserting new flights into the schedule; 
thus, a new four-engine flight is likely to be inserted into a schedule as a round trip between 
two points (and so makes a simple, but not necessarily the best, pairing), while two-engine 
flights may be inserted in a somewhat more circuitous fashion and so force new pairings. 

Thus, in the end, the part of the study that was of the greatest interest technically 
proved to have the lowest immediate payoff. Its greatest value lies in the fact that the airline 
now has a routine procedure that can be used in assigning flights to bases and so is protected 
against the promotion or loss of the very capable individual now doing the job. 


CONCLUSICN 

Probably the greatest over-all payoff from the study just described is the fact that it 
demonstrated to the policy-level personnel of the airline that existing empirical procedures 
may be routinized in a manner that makes possible an independent measure of the effectiveness 
ofthe procedure. A considerable portion of the final meeting, at which the report was pre- 
sented, was devoted to planning the next study—to be undertaken by the airline's own OR per- 
sonnel, Furthermore, some time was devoted to discussion of a much more complex study, 
to be undertaken in the not too distant future, which would probably have been considered 
impossible if it had been discussed before the results of the first study were in. 








ANNEX A 
The Assignment Problem Technique 


Lawrence Friedman and Arthur J. Yaspan 


In the text of this report, the statement is made that mathematical techniques do not 
exist for treating trip allocations for the airline system as a whole, but that a technique does 
exist for approximating optimal trip allocations for the whole system by treating it as a series 
of two-city problems. This technique is called the Assignment Problem Technique; it isa 
special case of the more general technique of linear programming. 

This annex describes the method of applying the assignment problem technique to the 
two-city problems into which the airlines system must be "broken" and "squared," as described 
inthe text. The first, second, and third matrices described in the text must be generated, either 
by hand or by a combination of hand and machine computation; then the third matrices may be 
"solved"' by the assignment problem technique. 

Two ways of solving the third matrices will be given in this annex. One is rigorous in 
the sense that it gives the optimum set of trip basings; this approach can be quite tedious and 
time-consuming, however, especially for large matrices. The other is an intuitive approach 
that requires more thought but less time on the part of the solver; it does not guarantee an opti- 
mum solution, but, properly done, comes very close to the optimum, 


FORMULATION OF THE ASSIGNMENT PROBLEM 

The assignment problem may be formulated as follows: 

Given n origins with a supply of 1 at each origin and n destinations with a demand of 1 
at each destination, with a cost, Cis of shipping 1 unit from origin i to destination j; then, if 
Xj equals the number of units to be shipped from origin i to destination j (5; = 1 or 0), the 
problem is to minimize the total shipping cost, C, where 


n n 
“ nia i ‘oo “ij “AY 


subject to 
(2) 


and 


(3) 
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Given a set of costs such as the one shown in Matrix 1, together with the requirement 
(Eq. 2 and 3) that the solution must contain one (and only one) element from each row and each 
column, then the optimum distribution of units is not immediately apparent. 
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METHOD OF REDUCED MATRICES (EXACT SOLUTION) 

Given a problem of this sort, the idea behind the method of reduced matrices is that, 
by simple and legitimate mathematical operations, the matrix may be reduced to an equivalent 
form in which the solution becomes obvious. It will be shown, for example, that Matrix 1 can 
be reduced to the equivalent matrix shown as Matrix 5 by steps that will be described. 






































Matrix 5 has the property that the optimal solution has a total cost of zero. If this 
matrix has a solution and if the total cost must be zero, then all allocations must occur in 
squares of zero cost. The allocations in Row 3 and Column 1 are immediately obvious, 
because each contains only one zero. The remaining allocations are easily determined by 
considering what remains to be allocated and the fact that the remaining allocations must take 
place in squares of zero cost. The solution is shown by diagonal lines in the upper left-hand 
corners of the squares entering into the solution. By referring back to Matrix 1, the total cost 
is seen to be 34, 

The procedure used in reducing Matrix 1 to Matrix 5 is as follows: 

Step 1: Take the smallest number in each row; subtract it from all numbers 
in its row. Matrix 1 then becomes Matrix 2. 
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Step 2: Take the smallest number in each column of Matrix 2. and subtract it 
from all numbers in its column. Matrix 2 becomes Matrix 3. 


: Matrix 3 consists of zeros and numbers. Draw horizontal and vertical 
lines so that lines pass through all zeros and so that the number of 
lines required to "cover" all zeros is minimum. Matrix 3 shows a set 
of lines that satisfies these requirements, For small matrices, the set 
of lines can be determined by inspection. A procedure for drawing these 
lines systematically for larger matrices is discussed below. 

: If the total number of lines required to cover the zeros is n, the size of 
the matrix, then the reduced matrix has a solution with cost zero and, 
as was shown earlier, the allocation becomes apparent. In Matrix 3, 
however, the total number of lines required to cover all zeros is three; 
therefore, a solution does not yet exist. 

: When a zero solution does not yet exist, examine the submatrix of 
squares that do not have a line through them (as shown by the submatrix 
enclosed by heavy lines in Matrix 3). Then select the smallest number 
in the submatrix, subtract it from each number in the submatrix, add it 
to all squares where two lines intersect, and leave all other numbers 
the same. By this procedure, Matrix 3 becomes Matrix 4. 


Step 6: Repeat Steps 3, 4, and 5 until a reduced matrix of zero cost is obtained. 
Matrix 5, discussed above, was obtained by this simple iterative pro- 
cedure and is a matrix of zero cost. 


SYSTEMATIC PROCEDURE FOR DRAWING LINES 

The procedure for drawing the minimum number of lines to cover all zeros is as fol- 
lows: 

First, after the first two steps in reducing the matrix have been accomplished, attempt 
a solution by marking as many zero cost squares as possible, subject to the restrictions. 
Matrix 3 is repeated below, with the attempted solution shown by the diagonal marks in the 
upper left hand corners, 
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Second, place an x at the foot of the columns with no marks and a check at the end of 
the rows with no marks, as shown above in Matrix 3. 

Third, in every checked row, search for squares containing zeros and check the cor- 
responding column, 

Fourth, in every checked column, search for zero cost squares that are marked and 
check the corresponding row. 

Fifth, repeat this procedure until either (1) a column marked with an x is checked or 
(2) no more checks can be made, If (1) occurs, the number of marked squares is not a mini- 
mum and improvements can be made by marking a square in that column and rearranging the 
marks, When an additional mark is obtained, the process must be repeated. Generally, one 
can determine the maximum number of marks on the first attempt. Where more than one 
possibility exists for marking squares, the choice may be made at random. If (2) occurs (that 


is, no more checks can be made), draw lines through each checked column and each unchecked 
row. 








Oo 








Matrix 3 





— 
wa 





3 
6 
3 
3 























x x 


This systematic procedure is usually unnecessary for small matrices, because the 
minimum number of lines is usually apparent. If a mistake is made in drawing the lines, the 
solution will not be affected; the number of iterations required for a solution may increase, 
however. 


THE APPROXIMATE METHOD (THE MAP-TACK METHOD) 

The intuitive approach to the solution of the assignment problem is sometimes called 
the map-tack method because of the convenience of using map tacks instead of pencil marks, 
inasmuch as map tacks can easily be moved from square to square as a solution is approached. 

In using the method, the procedure is to lay out the matrix on a sheet of paper and then 
to use map tacks to select a feasible solution. A good method for selecting a feasible solution 
is as follows: 

First, place a tack in the square of smallest cost. 

Second, place a tack in the square of smallest cost in the remaining matrix (that is, 
after eliminating the row and column in which the first tack was placed). 

Third, continue the procedure by placing the next tack in the square of smallest cost in 


the remaining matrix. This procedure might yield the feasible solution shown by the solid 
circles in Matrix 1A. 
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It will be noted that this solution has a cost of 36, as opposed to the previously determined 
optimum of 34. 

Fourth, examine pairs of squares containing tacks, searching for switches of two that 
will yield a lower total cost. A switch of two is defined as a change of two tacks in a manner 
that continues to satisfy the condition of one in each row and one in each column but results in 
a smaller total cost. Thus, in Matrix 1A, the dotted circles 4 and 10 have a smaller total than 
the corresponding circles 3 and 12. With this change, Matrix 1B (total cost 35) results, and it 
so happens that the additional switch of two indicated by the dotted circles reduces the total 
cost to the optimum (34). 








10 11 
6 oe 3 
Matrix 1B 13 12 7 

6 6 3) 
11 | 10] 12|@0| 8 


Fifth, in some cases (especially where the matrices are larger), one must search for 
switches of three, four, and more, until no more can be discovered. These higher-order 
switches are more difficult to find, and some experience is necessary before an effective 
performance can be achieved. 

Sixth, when no more switches can be discovered, the solution is called optimum, but it 
must be recognized that there may be a better solution that has not been found. Experience 
has shown that, with practice, matrices of the order of 20 x 20 can be solved in this manner 
with acceptable accuracy. 
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ANNEX B 
Calculation of the School Table 


Lawrence Friedman 


In the text of this report, a simple equation using what were called z values was devel- 
oped as a means of determining the size of classes in the Stewardess School. In this appendix, 
the means of arriving at the entries in the table are described. 

In the text, "surplus" and "shortage" were referred to; the method of computing opti- 
mum class sizes is based on the principle of minimizing the cost of these expected deviations 
from requirements, Here, the mathematical definition of these deviations is given: 

1, If a surplus is defined as the positive difference, at any given moment in time, 
between the number of available stewardesses and the number required for the system, and a 
shortage is the negative difference; and 

2. If the following symbols are used to mean: , 

z: Number of stewardesses available at the planning date, plus additions in the 
lead-time period, plus size of class being scheduled, minus requirements for 
the planning period; 

x: Losses during the lead-time period; 

y: Losses during the planning period; 

f(x): Probability that x losses occur during the lead-time period; 
g(y): Probability that y losses occur during the planning period. 
. Then, 
z-x (if positive) = surplus just after start of planning period; 
x-z (if positive) = shortage just after start of planning period; 
z-x-y (if positive) = surplus just before end of planning period; 
x+y-z (if positive) = shortage just before end of planning period. 

The objective, then, is to find that value of z that minimizes the cost of expected depar- 
tures from requirements for the planning period. The resulting optimum value of z can imme- 
diately be converted to optimum class size by following the instructions given in the text. 

In order to determine the z values, it is necessary to assume that the stewardess level 
during any given period (between class graduations) drops in linear fashion. Actually, as is 
shown in Figure B1, this loss occurs irregularly, in step fashion, but may be approximated by 
a straight line; there is little chance that significant error is introduced by this assumption of 
linearity. 
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Figure Bl - Assumed linearity of stewardess level in a period 


Inasmuch as the planner exercises discretion in determining the length of the planning 
period, it is convenient to express time in planning-period units. Thus, 


t=0 is the time at the start of a planning period; 
t=1 is the time at the end of a planning period; 
z-x-yt (if positive) is the surplus at time t; 
x + yt - z (if positive) is the shortage at time t. 


If z - x is positive, there will be a surplus at t = 0; this will persist until t = (z - x)/y or 
until t = 1, depending on which is smaller. If (z - x)/y < 1, there will be a shortage lasting 


from t = (z - x)/y until t < 1. Accordingly, the total departure-from-requirements cost will 
be: 
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The expected cost associated with the possibility 


z-x> 0, ~ee 


is then 





Cy (z - x)? a (x+y -2z)° 
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| f(x) g(y) dy dx. 
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Similarly, the expected cost associated with the possibility 


z-x>0, 1; == 


Z Z-X 2 
J j Cy x4 2 4 0 - xn - zy + xy f(g) g(y) dy dx. 


Finally, the expected cost associated with the remaining possibility 


z-x<s0 
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The total expected cost c(z) is the sum of these three double integrals. 
In order to obtain the value of z that minimizes c(z), we differentiate c(z) with 
respect to z and set this derivative equal to0. After some simplification, the result is 


Z 0 00 co 
(cy - &) i £(a)ax J __ ax)? SP ay + 09 i" £(x) ax i g(y) (22 - 2x - yay 


Z zZ-x 
+ (Cy - Cy) i ro) ax J g(y) (22 - 2x - y)dy=0. 


At this stage, it is not possible to solve explicitly for z unless C1 =Co; that is, unless 
the cost of a shortage per day is equal to the cost of a surplus per day. In the latter case, the 
first and third integrals drop out. Setting the second integral equal to 0 yields 


where X is the expected number of losses in the lead-time period, and y is the expected 
number of losses in the planning period. 


For the more realistic case where Cy 4Co, we rewrite the condition as 
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and inquire into the possibility of solving this equation numerically for z. 

To perform this calculation, one must make certain explicit assumptions about Co/c p 
f(x), and g(y). A likely value of 2 for Co/cy has already been advanced in the body of the text, 
As for the probability distributions of x and y, examination of resignation data by month 
indicated a significant seasonal trend and some yearly trend, but little evidence for any change 
in shape of the distribution once the mean was specified. It was decided to choose for each 
variable a rectangular distribution with range equal to 2/3 the mean; i.e., y equally likely 
from y = A to y = 2A with mean 3/2A, and x equally likely from x = B to x = 2B with mean 
3/2B. Formally, this amounts to setting 


f(x) = & (Bs xs 2B) 


f(x) = 0 (k<B or x > 2B) 
d 
g(y) = (As ys 2A) 
g(y)=0 (y <A or y> 2A). 
With these expressions for f(x) and g(y), and a value of 2 for Co/c}, condition * becomes 


1 1 re Z-xX ioe) d 
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where the upper limits may of course be truncated because of the assumed finite ranges of both 
the x andthe y distributions. (The numerical method developed for calculating z given 3/2A 
and 3/2B, is of no special interest.) 

The results of the calculations are the z-value entries in the table in the text; these 
may be converted to class size by following the procedure given in the text. 





CYCLIC SCHEDULING AND COMBINATORIAL TOPOLOGY: 
ASSIGNMENT AND ROUTING OF MOTIVE POWER TO MEET 
SCHEDULING AND MAINTENANCE REQUIREMENTS* 





A new method is developed for the nonlinear problem of optimal cyclic 
assignment of motive power, which derives essentially from the asso- 
ciation of a graph with the time-table. 











PART I 
A STATEMENT OF THE OPERATING PROBLEM 
OF THE FRISCO RAILROAD 


V. B. Gleaves 


St. Louis-San Francisco Railway Company 


Dieselization has been hailed as the greatest forward step taken by railroads in modern 
times. Diesel locomotives do not have to be cared for and pampered at each terminal as do 
steam locomotives. Except for requirements of fuel and water, they are capable of 5000 miles 
of continuous operation before coming into the shop for inspection. Apparently some railroads 
are not taking full advantage of the possibilities offered by the diesels. 

Prior to 1936 the average miles per day made by locomotives by all Class I railroads 
inthe United States was less than 100. This increased to 124.5 miles per day in 1943, prob- 
ably because of the war effort, then declined for several years, In 1952 there were more 
diesel units in operation than steam locomotives, and the average miles per day rose to 127. 

In 1955, with four times as many diesel units as steam locomotives, the average was only 147 
miles per day. 

Once a month Railway Age publishes statistics on the 51 largest railroads in the United 
Sates showing, among other things, average miles per locomotive per day. The February 18th 
issue shows that in November 1956 only the Lehigh Valley, KCS-L&A, Cotton Belt, Western 
Pacific, and Frisco averaged more than 200 miles per locomotive per day. Six of the railroads 
shown had less than 100 miles per day. The remainder averaged between 100 and 200 miles. 

In October 1956 the same five railroads were the only ones averaging more than 200 
miles per day. The Frisco was second highest, with 223.8 miles, exceeded only by the Lehigh 





‘Research for this paper was sponsored by the St. Louis-San Francisco Railway Company and 
Office of Naval Research 
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Valley with 252.2 miles. In September 1956 only three railroads exceeded 200 miles per day, 
counting the L&A and KCS as one railroad, and Frisco was second highest. 

Comparison of one railroad with another does not always give an accurate comparison, 
because the number of local trains, road switchers, and branch line trains has a great effect 
on the average miles per day. The Frisco Railroad, on the other hand, has its share of locals, 
switchers, and branch line service. It has 74 locomotives assigned to such service, with 42 
assigned to through freight service. If these locomotives were making only the average mile- 
age by U. S. railroads, it would require about 40 additional 3-unit locomotives; and a 3-unit 
locomotive costs half a million dollars. The railroad has been able to save this huge invest- 
ment by scheduling diesels on precision schedules and by having at Springfield, Missouri, 
around-the-clock transportation supervisors who handle diesel locomotives over the entire 
railroad. 

It would be pleasant to say that Frisco freight trains always operate on time and that 
locomotives never fail between inspection periods. But emergencies do occur, and the trans- 
portation supervisors make the necessary arrangements to provide locomotives. They also 
provide locomotives for the operation of all irregular extra trains, 

There are two "precision diesel pools" with 25 locomotives consisting of 93 units, 
protecting 11,114 train miles daily, which is an average of 445 miles per day. Not counting 
passenger units and switch engines, the Frisco owns 283 freight units. That means one-third 
of the total freight units are averaging 445 miles per day. Ninety units are assigned to local 
trains and road switchers and average only 102 miles per day. The remaining 100 units 
protect through trains that are not in precision schedules, irregular extra trains, and shopping 
margin. 

There are numerous terminals on the railroad that inspect switch engines, general 
purpose units, and passenger locomotives, but Springfield, Missouri, is the only inspection 
point for the road freight units on important through trains, In the precision diesel pools, each 
pool must begin and end at Springfield so that the locomotive can be cut out for at least twelve 


hours for inspection, after having made at least 5000 miles, and preferably not over 5500 miles. 


Each locomotive in a pool is scheduled to make this 5000 to 5500 miles on different trains all 
over the railroad and be gone from eleven to thirteen days before cutting out at Springfield. 
Since there are two precision pools, two locomotives are inspected daily at Springfield—one 
4-unit locomotive cutting out during the morning and another 4-unit locomotive in the evening. 
This permits having a uniform working force in the diesel shop. 

Some of the trains in these pools have 4-unit engines and some have 3-unit engines. 
The changes from 4 units to 3 units are accomplished by cutting one "B" unit out of a 4-unit 
locomotive, and when this is done the "B" unit must be scheduled to cut back in a particular 
engine so that it too will be back in Springfield at the regular inspection period. 

The Frisco has been scheduling diesel locomotives in this manner since the first road- 
freight diesels were purchased nine years ago. So you may wonder why the Frisco asked the 
Management Sciences Research Group at Purdue University to look into the problem of sched- 
uling diesels. Since there are twenty trains and twelve large terminals involved, it means that 
there are more than 100,000 different ways to make these diesel pools. Most of these 100,000 
combinations are unsatisfactory, because they require an additional locomotive, fail to bring 
an engine into Springfield for inspection, excessive mileage, insufficient mileage, close con- 
nections or other reasons, It is often necessary to try 100 or more combinations, requiring 
a week or more to find one workable schedule. But since there are probably a dozen or more 
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_ workable schedules, we do not know that we have the best possible schedule and do not take the 
time to do the months of work that would be required to explore all possibilities. 

Sometimes even a minor change in freight train schedules requires that diesel sched- 
ules be revised. This may happen two or three times per year. Sometimes a minor change in 
one freight train will not require a change in diesel schedules, but we find out several months 
later that this minor change permits a much better schedule. If it was not such a tedious task 
to revise schedules, we would be able to check our diesel schedules several times a year, or 
at least every time a minor change was made in any freight train. 

Another benefit we expect to receive from a fast, simplified method of making diesel 
schedules is a saving in locomotives assigned in local train service. With our present methods 
it would require several months to make a complete survey of the entire railroad. 

Quite often in a time-card meeting a proposal is made to change the time of a freight 
train. We need to know within a few hours if the proposed change will require additional, 
fewer, or the same number of locomotives, but that information is not now available on short 
notice. 

The problem, as it was presented to the Management Sciences Research Group, was to 
devise a simplified method of determining schedules for diesel locomotives that is fast enough 
to permit a clerk to determine several of the most workable schedules in a few hours' time so 
that the Transportation Department can pick the best schedule from the several possibilities. 
Each workable schedule must: 

(1) Require minimum possible number of locomotive units. 

(2) Provide for inspection of locomotives at Springfield, Missouri, at approximately 

5006 miles, but in no event to exceed 6000 miles. 


(3) Allow twelve hours for inspection, and stagger the inspection time of the two pools. 

(4) Avoid close connection for locomotives at terminals. 

(5) Avoid as much cutting in and cutting out of 'B" units as possible, and see that cut- 
out "B" units are provided a proper schedule. 








PART II 
GENERALIZATION AND ANALYSIS 


T. E. Bartlett and A. Charnes 


Purdue University 


I. INTRODUCTION 

In the initial stages of analysis a question was raised concerning the condition laid down 
by the Frisco Railroad that only the minimum number of units would be used to maintain the 
schedule. In a previous paper [1], an algorithm was developed for determining the minimum 
units required to maintain a fixed schedule of transport units. The expression which was 
arrived at (subject to certain constraints) appeared as follows: 


U=a+ DB; 
j ) 


minimum number of transport units 


total "runs" en route at beginning (or end) of 
schedule period 


maximum value of cumulative departures, less 
cumulative arrivals at terminal j during period. 


The algorithm did not present a method for allocating the units to the schedule, but 
merely determined the lower bound (without regard to feasibility) of the number of units to 
maintain the schedule. It is possible that for many schedules there exist a number of routings 
all of which require only the minimum number of units to maintain the system. We shall now 
proceed to develop a systematic method for obtaining these routings. It will be shown that 
some additional constraints which were not previously considered would serve to determine 
the choice from the available routings and thus provide an efficient procedure for generating 
feasible minimal allocation of transport units to maintain a fixed schedule. Also, a method 
will be developed for obtaining all routings which comply with restrictions on maintenance. 
These restrictions are given in terms of location of the servicing agency within the system, 
the amount of time necessary to service, and the required frequency of maintenance. 
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I, EQUIPMENT POOLS 

Current practice in a number of transport companies is to allocate equipment to 
"pools."" These refer to a fraction of the total runs to be accomplished during each period 
and consist of a fraction of the total number of transport units which are required to maintain 
the total schedule; the total of all pools comprise the total system. Within each pool the units 
are scheduled to accomplish in a specified sequence all of the runs assigned to the pool, per- 
forming each run in turn until all of the runs in the pool are accomplished. At the beginning 


and end of the sequence the transport unit is at the terminal at which a maintenance depot is 
located; the sequence having been so chosen that the time between arrival and departure is 
sufficient for maintenance. The maintenance requirement, in fact, largely determines the 
possible allocation of runs to pools. 

The maximum usage of equipment between maintenance is fairly rigidly determined by 
either legal requirements or operating characteristics. On the other hand, more frequent 
maintenance represents an increased cost with little payoff. The method of allocation by 
forming pools and distributing equipment to them offers an orderly procedure for both sched- 
uling the runs and ensuring that equipment will be maintained on a regular economic basis, 

An additional consideration which must be made in the formation of pools is the fact 
that certain runs call for equipment of a particular type. In other runs there may be a mini- 
mum specification for the transport unit which may be exceeded. In the concrete instance of 
the present study (dealing with railroads) this variation in transport units is limited to two 
types (either two "A" diesel locomotives, which have a cab and independent controls, coupled 
with one ''B" diesel locomotive, which has no cab and no independent controls, or two "A" 
locomotives coupled with two "B" locomotives). In practice, the runs are referred to as 4- 
unit and 3-unit runs, respectively. Ideally, pools are formed from runs which call for the 
same minimum equipment, but such is not always possible. Examination of a typical schedule 
will ordinarily show examples of possible cycles calling for low-cost equipment which do not 
connect with the maintenance terminal. This situation is resolved in scheduling of railroads 
by transferring the fourth unit among various 3-unit combinations. This procedure has a cost 
associated with it in terms of both time and labor dollars and is avoided when possible. 


Yet another consideration made in the formation of pools is the desirability of retaining © 


the same equipment on "through runs."" A "through run" is an-aggregation of runs which are 
regularly scheduled between two terminals and which traverse additional intermediate termi- 


nals en route. Any load which is transported on one run of the aggregate is usually transported 


over multiple legs of the combination; these aggregates are usually referred to as "trains" and 
as such are distinguished by unique numbers; the airline equivalent being "flights."" Arrivals 
and departures on these through runs at intermediate terminals are separated by relatively 
short time intervals. It becomes advantageous to assign the same transport units and all the 
runs in through runs whenever this is feasible, for the following reasons: 

(1) Layover at intermediate terminals is not usually sufficient to transfer equipment. 

(2) Service in the matter of freight deliveries may otherwise be greatly impaired. 
Since much of the freight to be transported over any leg of the through run usually is scheduled 
to go over the consecutive lags, a delay at any point will affect both the transport unit and the 
freight consist concerned in such a combination. Other assignments may result in a late unit 
bringing about a late departure for freight which is otherwise on time. 

In short, pools are made up so that as nearly as possible the following characteristics 
hold: 
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(1) Minimum capital investment in the form of transport units will be required. 

(2) The pools, in toto, will insure that all of the regularly scheduled runs are provided 
with suitable units and the scheduled arrival and departure times are not violated. 

(3) Sufficient time will be allowed between arrivals and departures, so that necessary 
transfers of single or multiple transport units may be made. 

(4) All transport units will be brought into the maintenance terminals after having 
been in use a large fraction of the maximum mileage permitted by legal or operating restric- 
tions. 

(5) The same units will be assigned to all the runs making up a through run or "train," 
whenever possible. 

(5) Transfer of the fourth unit is made only when absolutely essential. 

Examination of various real time-tables makes it obvious that it is impossible in many 
instances to provide all of the characteristics desired. It is necessary, therefore, to obtain 
some method for evaluating the violations in order to determine the actual assignment to be 
made. From discussions with operating personnel it was decided that the best procedure would 
be to hold the desired characteristics (1), (2), and (3) as firm constraints, place bounds on 
characteristic (4), and devise a systematic procedure for transcribing a number of alternate 
schedules within the remaining characteristics so that operating personnel may evaluate the 
entire schedule so provided. 


The development of the model for routing and allocation of transport units is thus based 
on the concept of pools. In order to proceed with the development it becomes necessary to 
examine the significance of layover or idle time of units at the various terminals in the system 


insofar as this affects the possible connection of runs with transport units. 

We shall ignore the possibilities of assigning transport units to trains at locations other 
than terminals, since such possibilities as a regular operating matter are usually not practical. 
Under this premise, it is basic to any method or routing and allocation to be able to specify 
and evaluate in an efficient and systematic manner the feasible transfers of units from run to 
run at each terminal. We shall find it convenient to develop such a system in terms of an 
"idle-equipment inventory." 


I, IDLE-EQUIPMENT INVENTORY 

At the time of arrival into a terminal, a transport unit is left idle at the terminal until 
itis removed by a departure. This unused equipment might be considered in idle-equipment 
inventory, which is increased by one unit for each arrival and decreased by one unit for each 
departure. It is, however, subject to various couplings of arrivals and departures which result 
in particular units being held idle for different periods of time. In the first paper of this 
Series, the inventory, in terms of equipment units multiplied by time units, was seen to have 
a lower bound, which was the sum of lower bounds at the various terminals that could be 
Separately determined for each from the time-table of runs. We shall now investigate the 
possible couplings of equipment consonant with utilizing the minimum number of transport 
units, As in the first paper of this series, it will suffice to consider three cases. The analysis 
of the three cases will in turn be reduced to the analysis of a single situation exhibiting the 
constraints essential to the problem at hand. The three cases will be presented in the same 
form in which they were stated in the previous paper. The graphical presentation, however, 
will be in terms of the inventory concept we have just introduced. 
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Incoming Runs: 7:00 AM,8:00 AM, 9:00 AM Case I: All departures in a typi- 

Outgoing Runs: 10:00 AM, 11:00AM, l:OOPM pical period at a termina] 
are later than all arrivals 
during this period. in 


Consider a schedule which calls ° 
for incoming runs at terminal 1 at 7:00 


a.m., 8:00 a.m., and 9:00 a.m. daily and 
outgoing runs from the same terminal 
leaving at 10:00 a.m., 11:00 a.m., and 
1:00 p.m, daily. If the horizontal axis 


gz. 3S 





Period K-1 Period K Period K+ 





Figure 1 - Idle-equipment inventory 


(example of Case I) refers to time and the vertical axis to 
idle-equipment inventory, representation 
may be made as in Figure 1. We shall adopt the notation which was used in the preceding 


paper: 8: 
th . ' :; ‘ ' beast 10 
ij - the i~ incoming run in chronological order during the period at terminal j , 

sj the jth outgoing run in chronological order during the period at terminal j 
tr 

ti the scheduled arrival time of ij 
& « A de 
ti the scheduled departure time of Tj ‘s 
tj - the scheduled arrival time of rij at the next terminal Ci 
m 
1 the length of the typical period in appropriate time units al 
the minimum total of idle times for transport units at terminal j " 
minimum total of idle times for transport units at all terminals tk 
minimum total running time for all transport units departing from terminal j r, 
el 
minimum total running time for all transport units departing from all terminals di 
total time of transport units required by a schedule Pl 
is 
minimum number of units to maintain schedule e 
r 
From the illustration we are able to infer clearly two significant points: 

(1) The total idle time for all units at the terminal in this case is invariant so long as ¢| 
equipment arriving during the period is utilized on (or "paired with") departures during the a 


same period. 

(2) Pairing any one or more of the departures with any of the arrivals from a previous 
period, or pairing any of the arrivals from one period with a departure during a subsequent 
period, only serves to increase the total idle-equipment inventory in terms of time multiplied 
by units on hand at a terminal. (The hatched area under the curve represents the total cost of 
idle equipment in time equivalent units.) 
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If all the terminals in a selected Incoming Runs: 7:00AM, 8:00AM, 10:00AM, 6:00 PM 
pool conformed to Case I, then every pos- CORES SO eee 
sible set of pairings of the first kind would 
imply utilization of the minimal number of 
transport units, 

Case II: Some departure is earlier 
than some arrival, but 
all arrivals and depar- 
tures within the period 
can be paired in at least 
one way so that each 
departure is later thanits 
paired arrival. 

Consider as an example of this case a schedule which calls for arrivals at 7:00 a.m., 

8:00 a.m., 10:00 a.m., and 6:00 p.m. and departures at 9:00 a.m., 1:00 p.m., 2:00 p.m., and 
10:00 p.m. In this case the schedule and idle-equipment inventory may be represented as in 
Figure 2. 

Again, as in Case I, we are able to infer some important relationships from the illus- 
tration. In the data at hand we may observe the following relationships: 

(1) The minimum inventory level is maintained if, and only if, the equipment for 
departing run 14 at 9:00 a.m. is obtained from either arrival Ti. at 7:00 a.m. or arrival To4 
at 8:00 a.m. A much more general statement may be made about this relationship. Specifi- 
cally, if we bear in mind the period 9:01 a.m. until 9:59 a.m. as an example of a "local mini- 
mum," it may be said that any departures leaving prior to a local minimum may be paired with 
arrivals which precede that departure. 

(2) The minimum inventory level in the case at hand is further restricted by the neces- 
sity to pair the departures To1> T34 at 1:00 p.m. and 2:00 p.m., respectively, with any two of 
the previous arrivals which are not paired with departure Tit 

(3) Finally, to maintain a minimum inventory level in the case at hand, departing run 
Ty, at 10:00 p.m. must be paired with arriving run r 41 2t 6:00 p.m. All these situations are 
encompassed in the generalization: To secure minimum idle-equipment inventory, * every 
departure lying between two successive local minima must be paired with an arrival which 
precedes the later minimum but is not prior to an antecedent minimum for which the inventory 
is zero, For example, if the earlier of the succeeding minima were at zero inventory, then 
each departure between these minima would have to be paired with an arrival between these 
minima. 

The local minima in any concrete situation are easily spotted, If the runs are sequenced 
chronologically, any departure followed by an arrival indicates a local minimum. For Case I 
and Case II the minimum is equal to zero whenever "Ky 1,j is the next later event after ee 
As may also be noted, in either Case I or Case II zero is always a minimum value and there 
are no negative minima, 





Figure 2 - Idle-equipment inventory 
(example of Case II) 





*Idle-equipment inventory will hereafter be referred to as "IEI."' 
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Incoming Runs: 7:OOAM, 1:00 PM, 3:00 PM 
Outgoing Runs: 8:00AM, 10:00AM, 4:00 PM 





Period K-! Period K Period K+! 





Figure 3 - Idle-equipment inventory 
(example of Case III) 


Incoming Runs: 7:00AM, 1:00 PM, 3:00 PM 
Outgoing Runs: 8:00AM, |0:00AM, 4:00PM 





Figure 4 - Idle-equipment inventory 
(example of Case III, modified) 


Case III: Up to some time during 
the period more depar- 
tures than arrivals have 
been scheduled. 

As an example of such a schedule, 
consider one which calls for arrivals at 
7:00 a.m., 1:00 p.m., and 3:00 p.m. and 
for departures at 8:00 a.m., 10:00 a.m., 
and 4:00 p.m. Our inventory plot then 
becomes that shown in Figure 3. 

To interpret the "negative inven- 
tory" of Case III, we must refer to the 
development in the previous paper where- 
in it is shown that it corresponds to a 
quantitative deficiency in connections 
possible within the period. In Case III, it 
became necessary both to "carry over"'a 
transport unit from a run arriving in the 
previous period and to carry over until 
the following period a transport unit from 
one of the incoming runs during the cur- 
rent period. This carry-over transforms 
the inventory as shown in Figure 4. 

The modified representation will 
have the same form as Case II, provided 
the beginning (and end) of the period 


shifted over suitably in time. Specifically, the new beginning may be chosen as any time at 
which the IEI is at its minimum. Such a time corresponds to the time of departure of that run 
which determines the maximum of cumulative departures less cumulative arrivals during the 
established time period. In other words, the times of minimum IEI are the times at which 8, 
of the preceding paper is attained. It is a useful convention to shift the beginning of the period 
to the first point in time during the original period when B is attained, or, what is the same 


thing, when the IEI is minimum, 


We are now in 2 position to enunciate a procedure for establishing in all cases the 
possible pairings which will not employ more than the minimum amount of equipment: 


1. Sequence chronologically the runs at the jt? 


terminal. 


2, Compute the cumulative excess of departures less arrivals during a single period; 
the maximum value so obtained is equal to B.. 
. Bp j is equal to zero, no shifting of the period is necessary. 
. £B j is greater than zero, shift the beginning of the period to the first point in the 
original period at which B j (minimum IEI) is attained. 
To secure minimum idle equipment inventory, every departure lying between two 
successive local minima must be paired with an arrival which precedes the later 
minimum but is not prior to an antecedent minimum for which the inventory is zero. 
The method above may be seen to embody all three cases, It will provide us with a 
simple procedure for obtaining from the information in time tables a set of restrictions which 
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will ensure that pools—and complete schedules—will utilize only the minimum number of trans- 
port units. The method so far is useful only when all transport units are of a standard type. In 
a later section this requirement also will be removed. 


Iv. FORMULATION AS A PROGRAMMING PROBLEM 

In the initial analysis, several mathematical programming formulations were explored. 
The form which appeared most promising to the concrete instance being studied (and which will 
be delineated in the following section) was a linear programming form in which the actual non- 
linear elements were to be borne in mind, but the hope was that the new type of constraints 
suggested by the preceding cyclical analysis would prove effective in yielding solutions which 
wouldalso satisfy the nonlinear conditions. 

In the initial formulation the problem was considered one of (a) assigning runs to pools 
so as to minimize the number of transport units while providing all runs with transport units, 
(b) insuring that the usage of each unit approximated the recommended mileage interval between 
maintenance, and (c) guaranteeing that each pool would have an idle period sufficient for routine 
maintenance. * 

Application of the "Minimum Transport Unit Algorithm" to some previous actual sched- 
ules, however, showed that they had, in fact, been employing minimum transport units. Also, 
by mathematical considerations (of "homology theory"), it was made plausible that one might 
expect very many possible minimal allocations to pools. Thus the emphasis was shifted from 
minimizing units to a criterion which might be more effective in reducing the possibilities for 
consideration, as well as in automatically satisfying the withheld nonlinear conditions. The 
objective selected was minimizing the sum of the absolute deviations from a criterion mileage 
related to the recommended maintenance interval. 

The notation adopted for the programming model would designate the lengths of the 
pools as say X, Y, Z; the "fraction" of the jth run at the jth terminal assigned to each pool 
being Xjjp Vij and Zi» respectively. The length of run Tij was aij: Thus the length of, or 
mileage traversed in, pool < would be: 


The length L of the entire system (which is generally different from the total trackage, 
because more than one run is scheduled over the same track) is the sum of the lengths of all 
of the pools, or 


(2) L=X+Y+Z. 


In order to write the functional for the system it was necessary to establish for each 
of the pools the mileage related to maintenance referred to above. It was a matter of policy 
that each pool was to discharge a unit for maintenance once during each period, regardless of 
the mileage in that pool. Since the possible losses which might be occasioned by breakdowns 
were considered to be roughly comparable for most parts of the system, a single figure was 
chosen for the criterion mileages of the various pools, e.g., L/3 for each of the three pools. 





*Additional constraints to connect ''through runs''and permit multiple-type transport units were 
considered but will not be discussed here. 
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Thus our objective would be to minimize |X - L/3| + |y-1L/3| + |Z - 1/3], 
subject to appropriate conditions. Recalling 


».4 


ij “4 
and the device of Charnes, Cooper, and Ferguson [2], this can be rephrased in linear form as: 
minimize x* +x +y*+y +z' +27, 


subject to © ayy Xj + xt -x"=L/3, 
i j 


etc. 


It was hoped that the nonlinear requirements that the Xia» Vip and Z; be either zero or 
that one might be assured by (1) using the additional linear requirements that the fractions of a 
run allocated to the different pools add up to one, e.g., 


Say * Faq * “59 * 


Xi + Vij + Zi =l1, 

(2) by procedures of the type employed in linear-programming computation of the travelling 
salesman's problem [4], and (3) by the following novel type of constraints which stem from the 
cyclic analysis of the minimum units paper and that analysis behind the IEI considerations. 

Rather than considering, as some [3] have done, a cyclic or periodic system to be 
adequately represented by a finite number of stages of repetitions of the period (when for 
equivalence an infinite number is required), the analysis of [1] and the foregoing has estab- 
lished, in consequence of cyclicity, certain critical points at which the events in a suitably 
transformed equivalent period may be laid out in sequential order. The new type of constraints 
are requirements intended to insure that only feasible connections preserving minimality of 
transport units are included in a pool. 

To illustrate, such constraints for the example of Case I above would be (renaming 
x's, y's, and z's for outgoing runs as x's, y's, and z's): 


a 


X41 + XQy + Xgq - Xyq - Xqy - Xgy = 0 


¥41 + Yoy + ¥31 - 41 ~ You - ¥g1 =9- 
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For the example above of Case II, the constraints (with similar renaming) would be: 
X11 ~ X41 - XQ <0 


X91 + X3y ~ Xyy ~ XQq ~ Xgq < 9 


X4y - Xq, = 0- 


Constraints of a similar character covered the example in Case III. 

The special character of the maintenance terminal, in requiring a layover of sufficient 
time to insure maintenance and in requiring the apportionment of the possible "layover" pair- 
ings to the various, was reflected by conditions such as: 


Finally, non-negativity constraints were placed on all of the variables. 

In the prototype model twenty runs were to be allocated to two pools. 

For the real problem no standard linear programming calculations were ever made on 
the problem exhibited, for it was noticed almost at once that it was possible to trace out a 
feasible closed path (in space and time) by starting in a row of the matrix on a plus 1, going to 
a - 1 in the new row, and so forth, back to a row referring to a condition at the starting termi- 
nal. The question then posed was what does this section of the programming matrix represent ‘ 
Our answer, which we shall partially elaborate in the succeeding sections and more fully 
present elsewhere, lay in establishing its connection with combinatorial topology. The sub- 
matrix is a type of simultaneous presentation of a collection of the "incidence" matrices of 
certain "chains" and nodes in "'graphs,"' which we can associate with every system of the type 
we have been considering. 

Alternatively, one can associate with each system a single graph from which links are 
to be struck off to form the "pool"’ subgraphs of the system, (Since this work was done, a paper 
by O. Ore has appeared [5] which may have application to the question of existence of "'pools."') 


V. THE SCHEDULE AS A GRAPH 

By an oriented or directed graph is meant a collection of "nodes" plus "links" connect- 
ing (or incident on) nodes with assigned orientation or direction. In our situation a convenient 
representation is: 

Node =a specific time at a specific terminal. The time is that of a 
local maximum of the IEI, 

Link incident at a node = arriving run prior to a node time (or departing run after a node 
time) which may be paired with a departing run after a node 
time (respectively, arriving run prior to a node time) without 
violating minimality of equipment. Such time portions of 
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arriving runs are links positively incident on the node and 
such time portions of departing runs are links negatively 
incident on the node. 

The total of an arriving or departing runs is thus a "chain" of links. 

We set up the "incidence matrix" of the nodes of our system with ''chains" or links 
corresponding to runs. To each node corresponds a row of the matrix and to each run cor- 
responds a column of the matrix. A plus one (or, respectively, a minus one) is entered if the 
chain has a link positively (or, respectively, negatively) incident on the node. Zeros or blanks 
are inserted where there is no incidence of chain and node. For example, we consider some 
actual time-table data for a portion of the system under consideration. 


Time-Table for Prototype Problem 
136 835 836 


1430 2100 lv StL arr 700 
600 630 arr Mem lv~- 2000 
330 
1715 
1430 
430 


Total mileage scheduled: 4816 (out of 10,714 actual) 


The first step necessary in setting up the incidence matrix is as follows: 
(a) Arrange the arrivals and departures at each terminal in chronological order. 
(b) Assign a minus one to the departures and a plus one to arrivals. 
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(c) Compute the cumulative total of the assigned values in their chronological sequence 
as follows: 


Determining the 8 for Kansas City Terminal 
Prototype Problem 


Time of Day Event Value Cumulative Total 





1000 Dep -1 -1 
1159 Arr +1 0 
1430 Arr +1 +1 
2020 Dep -1 0 
2140 Dep -1 -1 
2359 Arr +1 0 


This computation gives us two things: It gives the value of 8, for determining the 
minimum number of units to maintain the schedule, and it also provides the information for 
the transformation of the time period so that the nodes of our graph will all fall later in the 
sequence than their associated positively incident links. 

In step two of the computation, the arrivals and departures which precede the first 
event after that which determines the smallest cumulative total are shifted to the end of the 
sequence in the same order. In the example for Kansas City, there are two times at which the 
smallest cumulative total occurs: at 1000 and at 2140. In such cases, either point may be 
selected as the basis for the shift, although the first is here arbitrarily selected. After the 
transformation, step one above is repeated as follows: 


Transformed Sequence of Arrivals and Departures 
at Kansas City Prototype Problem 


Time of Day Event Value Cumulative Total 





1159 Arr +1 +1 
1430 Arr +1 +2 
2020 Dep -1 +1 
2140 Dep -1 0 
2359 Arr +1 +1 
1000 Dep -1 0 


The next step is to determine the nodes of the matrix and the runs which are negatively 
and positively incident on the node. In order to accomplish this in a simple manner, place a 
Single line under each arrival which is followed by a departure and a double line at the beginning 
of the series under each departure for which the cumulative total is zero. Each single line then 
represents a node. All of the runs which are incident to a particular node lie between the double 
lines, All such arrivals prior to a node or single line are positively incident on the node, and 
all such departure following the node or single line are negatively incident on the node. 
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The final step in the development of the incidence matrix is to define the rows and 
columns of the array and make the appropriate entries. Each of the columns corresponds to 
one of the runs in the system, and each of the rows corresponds to one of the nodes, there 
being at least one node for each terminal, For the row which represents a particular node, a 
plus one (or minus one) is centered under each column if that run is an arrival positively 
incident (or, respectively, a departure negatively incident) on the node. Figure 5 represents 
the incidence matrix for the prototype problem. 

In the matrix as it now stands it is possible to determine a "pool" in a very simple 
manner. Connect a minus one with a plus one in the same column, connect that plus one with 
a minus one in its row, and repeat, connecting plus ones and minus ones alternatively within 
the same column and then the same row. Only one pairing or connection of entries within a 
column is permitted. Each entry, if used, is connected with precisely two other entries, one 
in its row and one in its column. When the first or starting entry is connected the second 
time, a cycle or pool is formed which may be interpreted as follows: If a locomotive is 
scheduled to make the run which corresponds to the first entry, a feasible next run for the 
locomotive is the next column connected with the first run, and so on. The pairings of runs, 
furthermore, preserve minimality of equipment. 

If the pairings which are chosen form two cycles, each of which have ten runs in the 
prototype problem, they will each be sufficiently close the desired mileage between mainte- 
nance, The shortest ten of the twenty runs in the prototype total 2086 miles, or 317 miles 
short of one-half of the total mileage. Actually, it is impossible to make a cycle of the ten 
shortest runs, The fact that terminals in the Frisco railroad are approximately the same 
distance from each other is not unusual in transportation systems. Other operating require- 
ments usually make this a characteristic of the system. It is ordinarily necessary, therefore, 
to have included in a pool some fraction of the total number of runs to insure that the proper 
mileage be included in a pool. In the larger problem for the Frisco there are fifty runs with a 
total of 11,114 train miles daily, or an average of 222 miles each, If a cycle consists of twenty- 
five legs, there is an extremely high probability that a pool randomly formed will have less than 
6000 miles and more than 5000 miles. 

In order that each pool have a layover in the maintenance terminal, it is only necessary 
that some pair of departing and incoming runs at the maintenance terminal which provides such 
time be included in each pool. These connections in the real problem are those which provide 
a twelve-hour layover at the Springfield terminal. A simple method for providing the proper 
maintenance time is to select arbitrarily two such pairs of arrivals and departures and assign 
one pair to the first pool and the other pair to the second pool. 

"Close connections" may be avoided in general by the following procedure: 

(1) Adding a constant to the arrival times at terminals sufficient to cover the expected 
operating contingencies. This is done prior to computing the chronological sequence at each 
terminal, 

(2) Arbitrarily excluding one of a pair of oppositely incident runs which are too close 
in time as incident on a particular node. If this is done, it is necessary to set up an additional 
node incident with the run excluded and those runs which it is practical to connect with it. 

Finally, in order that cutting in and out of runs be avoided, it is possible to utilize a 
simple device* which consists of defining the entrees as "plus four," "plus three," "minus four," 





*This device permits, among other things, the application of the technique to systems in which 
a 4-unit locomotive may be used also as a pair of 2-unit locomotives. 
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and "minus three," rather than plus one and minus one. The values used indicate the number 
of units required on the run. In order to reduce the number of times a fourth unit is cut out, 

it is only necessary to pair plus fours with minus fours and plus threes with minus threes 
whenever possible. When it is not possible to pair a three with a three, a pairing of a plus or 
minus four with a minus or plus three indicates that a unit is cut in or out. While the analysis 
of this device is not complete, it appears from the data at hand that it will lead to the minimum 
number of cut-outs and cut-ins for the system. 


VI. CONCLUSIONS 

The use of the incidence matrix and the related computational devices indicated permits 
us to solve essentially the problem which was presented by the Frisco Raikroad. As is to be 
expected from such a model, the analysis has uncovered many interesting features which were 
not obvious at the start. Additional analysis is currently under way, which should lead to an 
even more useful model that may have much wider application, 
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ABSTRACT 


The many similarities between communication and transportation network problems 
make it apparent that a closer association is desirable between analysts in the two fields. In 
both fields the problems involve the adequate description of the network itself (capacities, etc.), 
the demands of the users (traffic intensities, etc.), and various economic factors. Two general 
classes of problems are readily discernible: problems of most efficient utilization of an exist- 
ing network and problems of optimal plant extension, both of which involve large numbers of 
highly interconnected variables. 

This paper deals with techniques, useful in treating communication network problems, 
which may have similar utility in the transportation field. 

1. Dating from the time of Erlang [5], probability theory, in particular queueing 
theory, has been extensively applied to traffic analyses, routing problems, and equipment 
design problems [13,14,16]. 

2. To deal with those network problems that are analytically intractable, use has bee. 
made of simulation devices known as "throwdown machines" [6]. Present-day advances in 
large-scale high-speed. electronic computers hold even greater promise for future extensions 
of this technique. 

3. Certain network problems at the system level are amenable to treatment through 
the use of linear programming. These include various routing and design problems [8,9]. 

4. Graph theory [11], a branch of topology, is useful in characterizing networks pei se 
and in dealing with certain traffic problems [15]. A result of Kruskal [12] can be used in an 
optimal design problem. 

5. Boolean algebra is useful in the description, manipulation, and design of two- 
terminal switching networks [17,10]. 

6. Dynamic programming [1,2] holds promise for applications in communication 
system problems [9]. Of interest to transportation researchers is the fact that Bellmai has 
applied it to certain Hitchcock-Koopmans transportation problems and to a railroad routing 
problem [3]. Boldyreff's flooding technique has been used in this connection also [4]. 

Progress made by communications analysts along other lines, particularly in the use 
of computers for control as well as operations research, is indicated. 
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Analysis of cyclic-linked cargo-handling systems is presented and 
solved by Monte Carlo techniques for several typical instances. 











INTRODUCTION 

The Maritime Industry is a highly complex transportation system. It embraces a 
diversity of interests and equipment, a wide variety of skilled and unskilled manpower, and 
a collection of heterogeneous facilities. The loading and unloading of ships, commonly 
referred to as cargo handling, is the heart of this system. So much time is required for 
cargo handling that it is a key factor in determining the number of voyages a vessel can make 
per year. 

The design of methods of handling cargo is largely an art. It will not become a science 
until the unformulated experiences are assembled and the physical laws which govern the 
movement and handling of the individual packages are deduced. While water flowing through 
a pipe and a ball rolling down an inclined plane are useful concepts, such analogies are not 
adequate for quantitative calculations of performance. There is, however, a widely accepted 
measure of merit for cargo-handling systems. It is the rate of flow expressed in either weight 
or volume tons per hour. The average rate of flow, commonly called "productivity,"' is used as 
a means of comparing different methods of work, different facilities, etc. Productivity is also 
the basis of planning operations, and is the key to all costs. 

A limited amount of productivity data is available. It ranges from broad over-all 
averages to figures for loading specific case goods commodities. For most commodities, 
productivity data are few, and those which are available are variable and controversial. Hence, 
interpolation and extrapolation are hazardous. To facilitate the design of new cargo-handling 
systems and the evaluation of existing ones, over a wide range of parameters, a rational tech- 
nique of prediction based upon a theory of flow for discrete entities is needed. 

This paper presents the analysis of productivity that has been predicated upon the 
idealizations of the materials-handling system developed by the Staff of the Cargo-Handling 
Research Project. Research, sponsored by the Office of Naval Research, has been in progress 





*Research on this paper was sponsored by the Office of Naval Research. 


223 





224 R. R. O'NEILL 


on the Los Angeles campus of the University of California since 1951. Studies have been 
directed toward the determination of relationships between the basic elements of the gen- 
eralized cargo-handling system, with the requirement that these relationships should provide 
criteria for both the design and evaluation of specific systems. Furthermore, the program 
has had as its objective a study of the entire system from shipper to receiver. Although 
analytical studies must be made on the individual components or subcomponents of any sys- 
tem, it is not until the results of analysis are used to synthesize an effective operating system 
that tangible improvements are realized. 

The study of productivity had as its point of departure observations of actual loading 
and unloading operations in several ports in the United States. It was observed that the loading 
and unloading of cargo involves a great variety of shapes and sizes of commodities. Actually, 
cargo is a composite of commodities from many warehouses and factories all over the world. 
Hundreds of shippers and consignees are represented in the commodities loaded and discharged 
at the several ports of call of a ship. The loading and unloading takes place on pier facilities 
which vary widely. Some are antiquated, while others are as modern as any industrial struc- 
ture. Furthermore, there is no standard design for a ship. Finally, the methods for handling 
cargo are not the same from port to port; the number of men in the "gang,"' which unloads a 
single hatch, varies from port to port, because gang size is negotiated between the employer 
and the union. In Los Angeles the standard gang is eighteen men, in New Orleans nineteen 
men, and in New York City twenty-three men. 

Because cargo handling varies so, it is not possible to describe a standard loading and 
unloading operation, but the following description will illustrate the basic operations. In Fig- 
ure 1, cargo is taken from its position of final rest on the pier and transported by the fork-lift 
truck to the apron opposite the hold of the vessel into which it is to be loaded. After each load 
is delivered to the apron, the fork-lift truck returns to the transit shed for another load. The 
cargo is picked up from the apron by the hook, hoisted up and over the side of the vessel, 
lowered into the hold, and detached from the hook. The hook then returns to the apron for 
another draft. In the hold each package is taken from the pallet by one or more hold men, 
depending upon the size of the unit being handled, and stowed. While methods differ from port 


Figure 1 - Schematic diagram of loading maritime cargo 
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to port and commodity to commodity, general cargo is moved piece by piece from the pier to 
the ship, or vice versa, by a sequence of movements and handlings that are repeated at regu- 
lar intervals. 


IDEALIZATION OF CARGO HANDLING SYSTEM 

The handling of cargo has been "idealized,""! as shown in Figure 2. The factors of the 
system are Classified into a limited group of independent elements: (1) the facility, (2) the 
inland transport, (3) the sea transport, (4) the commodity, and (5) the control. The sixth ele- 
ment, the process, is dependent upon the others and is the path of the transporting agent and 


the commodity. 
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Figure 2 - (a) The cargo-handling system, (b) cargo 
handling as a series of cyclic operations 


As the Figure 2 suggests, the process of loading and unloading of cargo may be con- 
sidered as a sequence of cyclic operations. In each segment of the operation there is a trans- 
porting agent which is a self-contained unit capable of performing all of the operations neces- 
sary to transport the commodity from one place to another. The round trip path of this 
transporting agent is considered to be a "link"; one round trip is a "cycle."" Note that there 
may be more than one transporting agent associated with a particular link; for example, there 
may be four men in the hold. The region where the cargo is transferred from one transporting 
agent to another is called a “anode.” Thetransfer may or may not be direct, and the commodity 
may remain at a node for any interval of time. 

If during the loading or unloading operation there is no net storage at the interior nodes, 
it can be readily seen that the productivity of the entire system is the same as the productivity 
of any single component. Thus, attention may be directed toward the productivity of the indi- 
vidual links, In each link the transporting agent picks up the load, moves it to the next location, 
releases the load, and returns for another. While engaged in these productive activities, how- 
ever, the transporting agent may encounter delays. There may be nothing to pick up; there 
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may be a stoppage during an actual operation or move- 
ment; or there may be no place to put the load. All of 
ieee these delays can be grouped into two categories: "induced" 
DELAY and "internal."' Induced delays arise when the transporting 
’ gga agent in the adjacent link either does not provide a load for 
pickup or is not available to receive the load when it 
arrives at the node. Internal delays are all stoppages, 
avoidable or unavoidable, during the actual operation or 
movement. The total time that the transporting agent 
TOTAL TIME takes to complete one round trip can be represented as 
shown in Figure 3. 
rs tr ecm pedbaad It is desirable to calculate the productivity in terms 
of the working times alone. At first it may appear that it is 
necessary to examine only the slowest link. Here, there 
should be no delay other than internal. Field studies [2], however, have shown that even the 
slowest link encounters considerable induced delay. This is the result of the cycle-to-cycle 
variation in the time required to perform each basic operation. If the transporting agent in the 
slowest link makes a fast trip at the same time that the transporting agent in one of the adjacent 


links makes a slow trip, the slowest transporting agent will have to wait, and thus encounter an 
induced delay. 


The average time per cycle for several 
transporting agents operating in sequence is 


DELAY CAUSED BY shown graphically in Figure 4, The time bars 
VARIATION are divided into three segments: (1) working 
time (2) delay time caused by the work methods 
and imbalance among the capacities of the trans- 
porting agents, and (3) delay time caused by 
cycle-to-cycle variation. Because there is no 
net storage at the nodes, the average time per 
cycle must be the same for each link. But the 
average working time for each transporting 
agent is not necessarily the same. Calculation 
of the delay caused by unbalance is relatively 
simple, provided the necessary methods study 
data are available. Calculation of the delay 
7 oe & caused by cycle-to-cycle variation, unfortu- 
LINK nately, is not so simple. 
? f The magnitude of the cycle-to-cycle 
a variation is shown by Figure 5. It is seen that 
the time required to perform each basic oper- 
ation is subject to unknown factors and may 
have different values, This introduces the notion of stochastic variation rather than deter- 
ministic behavior and the concept of average or mean delay. 
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CALCULATION OF PRODUCTIVITY 
Equations have been written which take into account the cycle-to-cycle variation, and 
for the simpler cases they have been solved [3]. The results for the simplest case are shown 
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Figure 5 - Example field data for hook link 


in Figure 6. The system consists of two links, with one transporting agent in each link. One 
unit of commodity is handled by each transporting agent. This unit is transferred directly 
from the first transporting agent to the second. Furthermore, the working activities are 
combined so that there are only two: (1) pickup plus transport loaded and (2) release plus 
return empty. It was possible to express the induced delay in this system as an explicit 
function of the working activity times; hence the average delay could be calculated for par- 
ticular density functions. The induced delay is shown as a function of the unbalance between 
the two links for two distributions of element working times. It is interesting to note at this 
point that the delay increases as the two-link system approaches balance. 

The three-link, zero storage system, shown in Figure 7, has also been studied ana- 
lytically. Although recurrence equations could be written, the delays couldn't be expressed 
explicitly in terms of working times. Consequently, it was not possible to calculate the mean 
delay, although upper and lower bounds were established, as shown. 

An impasse had been reached. It was necessary to use another technique, simulation, 
to obtain quantitative estimates of mean delay time for more complex and realistic systems. 
Simulation, however, involves point-by-point computations. When there are many parameters, 
the number of points required to establish the functional relationships between all of the 
parameters is very large. After careful consideration, it was decided to carry out the com- 
putations in such a manner that the relative influence of each individual parameter could be 
Studied. One parameter was varied at a time. The parameters involved in the link-node 
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Figure 6 - Average induced delay for two-link system 


idealization are six in number: (1) the number of links in the system; (2) the arrangement of 
the links, that is, whether the links are arranged in series, parallel, or some combination 
thereof; (3) the number of transporting agents operating in each link; (4) the relative size of 
the unit of commodity handled by each transporting agent per trip; (5) the number of units of 
commodity that may be stored temporarily at the nodes; and (6) the distributions of the pro- 
ductive activity times, P (pick-up), L (transport loaded), R (release), and E (return empty). 
Simulation was accomplished by adding together in proper sequence the particular 
values of the working times which were selected at random, allowing for the induced delay 
time as it occurred. In this way, the time required to move a single unit of commodity from 
its initial point of rest to its final point of rest was calculated. By repeating these calculations 
a large number of times, a good estimate of the average productivity, including delay time, was 
determined, The calculations were carried out on SWAC, the high-speed digital computer at 
the Numerical Analysis Research Project on the University of California, Los Angeles, campus. 
The simulation computations themselves may best be illustrated by an example. Con- 
sider Figure 8. The movement of material from A to B is to be accomplished in three steps. 
There are three links with direct transfer between links. There is one transporting agent in 
each link; it moves one unit of commodity at a time. At the start all transporting agents are 
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Figure 7 - Average induced delay for three-link system 


ready to beginthe return emptyphase. Simulation involves 
the random selection of element times. This was accom- 
plished by storing the distributions of element working 
times in the high-speed memory of SWAC, Recall from 
Figure 5 that the distributions of element times are dif- 
ferent. In each case, however, there is some minimum 
time required to perform the operation or movement. 
Furthermore, each distribution has a single modal value— 
one time that occurs most frequently. The distributions 
are skewed to the right because they include the internal 
delays in the right-hand tails. Rather than by the use of 
one particular set of observed element times, the calcu- 
lations were based upon synthetic distribution functions 
which had, approximately, the same mean values and 
Standard deviations as those observed in the field. The Figure 8 - Simulation of flow 
distribution functions are shown in Figure 9. The shapes through a three-link system with 
of two of the synthetic distribution functions do resemble Grom tremmer tetwess ane 


(one transporting agent and one 
those of the field data, but, more important, they may be unit of commodity in each link) 
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UNIFORM considered as probable limiting cases within which 
the actual distributions lie. The use of the expo- 
nential density is equivalent to supposing that the 
activity time is the sum of the absolute minimum 
time, corresponding to the most favorable circum- 
stances, and an additional random component of 
EXPONENTIAL time. However, the modal value of the exponential 
density is the minimum value. Since this may seem 
unlikely to some, the log-normal distribution was 
used for comparison. The uniform density was 
taken as an extreme case. 

A sequence of random numbers was gen- 
erated by SWAC through the use of an arithmetical 
process. These random numbers were used to 
select particular values of the prestored activity 
times. Values of E, P, and L thus selected are 

Figure 9 - Comparison of three added successively to the cumulative time for the 
a transporting agent in the first link, Ty. Because 
no storage is allowable, the transporting agent in 
the first link cannot release the first unit of com- 
modity until the transporting agent in the second link is in position. Therefore, another value 
of E is selected at random and added to the cumulative time for the transporting agent in the 
second link, To. At this point, Ty and Ts are compared and both are made equal to the 
larger. Next, another value of P is selected. It is added to both T, and To, because the 
pickup time for the second transporting agent is considered to be equal to the release time of 
the first. L is selected at random and added to To. E is selected at random and added to To. 
Then Ts and Ts are compared. This process is continued until the first unit of commodity 
is moved through all three links. At the end of the first round of calculations Ty, To, and T, 
represent the cumulative time after the transporting agents in the first, second, and third links 
have successively released the first unit of commodity. This process was repeated over and 
over. Ts therefore increases by increments as the first, second, third, etc., units of com- 
modity arrive at B. The basis of this method is sufficient repetition so that the average of the 
individual values becomes an acceptable estimate of the true average. In order to establish 
the number of times that the calculation should be repeated, a pilot run was made, and it was 
concluded that the average, based on 2500 units of commodity, would be satisfactory. 
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RESULTS 

In this manner a number of different systems were studied, making it possible to 
describe the flow of material as a function of the six parameters which determine the pro- 
ductivity of the system [4]. The results have been expressed in terms of "effectiveness," 
which is the ratio of the productivity of the entire system to the productivity of the individual 
components. The effectiveness, e, was calculated as the ratio of the time required for n units 
of commodity to move through the system when there are no induced delays to the time which 
was calculated, allowing for delays, T(n). 
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What appear to be the more significant results follow: 

Inasmuch as the movement of cargo has been considered as a stochastic process, it is 
to be expected that the density or distribution functions of the element working times are fun- 
damental quantities. Recall that three different density functions, exponential, log-normal, 
and uniform, were used, Their influence may be observed by plotting a family of curves for 
a particular system with the coefficient of variation (ratio of standard deviation to the mean) as 
the abscissa and the effectiveness as the ordinate. The calculations for systems of two and 
then twenty links are shown in Figure 10. The striking observation is the relatively minor 
influence of the shape of the distributions, as characterized by the third, fourth, and higher 
moments, 

Values of the coefficient of variation ranging 
from zero to unity were used in the calculations. 
There is some evidence to indicate that the coefficient 
of variation is approximately one quarter for the 
operations involved in cargo handling, but the calcu- 
lations were extended in order to investigate the 
influence of variability. Furthermore, there are 
many materials-handling operations which are more 
variable. Figure 11 shows some observed values of , : eenene. ceuerry 
standard deviation and average cycle time for cargo o eee eer 
handling. 

As expected, the number of links comprising 

O01 02 03 04 05 06 07 O08 O9 10 
the system is found to be a significant parameter. COEFFICIENT OF VARIATION 
Induced delay is caused by lack of synchronization. 

As the number of links is increased, the probability A806 1 Indsuence oF distribu 
of asynchronization between two adjacent links also no storage) 

increases. Thus, the average induced delay is 

increased, and the effectiveness of the system is 

decreased. The relationship between the number of links and the effectiveness was determined 
for a balanced system of N-links connected in series, balance referring to the ability of the 
transporting agents to move the commodities at equal rates. In each link a single transporting 
agent carried a single unit of commodity which was transferred directly from one transporting 
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Figure 11 - Relationship between mean and standard 
deviation for one group of field observations 
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Figure 12 - Influence of number of links 
(no storage, one transporting agent) 


agent to the next. The results are shown in Figure 12, The general shape of the curve may be 
explained intuitively as follows: Whenever transporting agents in adjacent links have fast, 
average, or slow cycles simultaneously, the induced delay is zero or almost zero. A delay 
will result only when one transporting agent completes its cycle in less than average time and 
must wait for the other transporting agent that has taken longer to complete its cycle. But 
when several transporting agents are delayed simultaneously, the total delay is no greater than 
when only the first transporting agent in the sequence is delayed, It, therefore, seems reason- 
able that the induced delay is caused principally by the transporting agents in the first few 
links, accounting for the general shape of the curve. For systems of few links the effectiveness 
can be increased by reducing both the number of links and the variability in the time required 
to perform the individual operations and movements, For systems of many links (e.g., bucket 
brigade) only changes in the cycle-to-cycle variability will influence the effectiveness. 

The single factor which contributed to the induced delay of the transporting agents in 
the simple N-link system was the requirement that each unit of commodity must be transferred 
directly from one transporting agent to the next. In other words, both transporting agents had 
to be at the node before either could perform the release or pickup operation. This delay would 
be reduced considerably if either the transporting agent arriving with the load could set it down 
without delay and continue on its cycle, or if the transporting agent arriving empty could pick 
up a unit of commodity from a stock pile. In the limit, this delay could be eliminated com- 
pletely by providing a very large stock pile at each node. For this phase of the study a three- 
link system was selected. There is no difficulty in accounting for the general shape of the 
curves shown in Figure 13, from intuitive considerations, The value of effectiveness would be 
expected to have a minimum for conditions of zero allowable storage and increase monotoni- 
cally to an asymptotic value of unity as the number of units of allowable storage increased. It 
is significant to note that a substantial increase in productivity can be achieved by making 
provisions for only one or two units of temporary storage. 

The single cause of delay that has been considered up to this point is the process of 
transferring the commodity at the node. If, however, more than one transporting agent is 
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Figure 13 - Influence of amount of storage 
(three links, one transporting agent) 


operating within a single link, a delay may be encountered, because the transporting agents 
must maintain their positions with respect to each other. For example, in the hold of the ship 
two or more men may perform the same task of carrying one carton from the pallet being 
unloaded under the hook to the wings where it is finally stowed. In this case a queue may be 
formed either at the pallet or at the place where the cargo is being secured. Again a three- 
link system was selected for study. The results are shown in Figure 14. The curve for the 
coefficient of variation equal to zero breaks at M = 4, because the first transporting agent 
delivers one unit of commodity and returns for the next pickup in less time than that required 
for the other transporting agents to make their first pickup. This is a consequence of the 
particular values assumed for P, L, and E. In the terminology of queueing theory, M = 5 
corresponds to the condition where the utilization factor* is greater than unity. In this case, 
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Figure 14 - Influence of number of transporting agents 
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*The utilization factor is the ratio of the average service time to the average interval between 
arrivals. 
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the results are sensitive to the proportion of the total cycle time that is required for pickup 
and release, In any event, the interaction among transporting agents within the same link is 
another factor to consider when deciding whether to use one large piece of materials-handling 
equipment or several smaller units. 

By varying one parameter at a time and using relatively simple balanced systems, the 
influences of the individual parameters have been examined. There are many actual situations, 
however, where the arrangement of the links is more complex and where the transporting agents 
are not matched. In order to demonstrate the power of simulation applied to the flow or transfer 
problem, the calculations were carried out for a more typical cargo-loading operation. The 
methods and time data were observed on the Los Angeles waterfront [5] during the summer of 
1955. The generalized activity chart, shown in Figure 15, indicates simultaneous activities 
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Figure 15 - Generalized activity chart, maritime cargo-loading 
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horizontally and sequence in time verti- 

cally. The single fork-lift truck is the 

transporting agent on the pier. It moves 

two loaded pallets per trip to a spot onthe 
apron under the hook. The hook lifts the Ric. 
pallets one at a time, servicing the two J ALTERNATE CvcLES 
hold gangs withalternate cycles, The first 

pallet is taken to a team of two men work- 

ing in the starboard wing of the hold. The 

loaded pallet is left there and the hook is 

then swung over to the port side where it 

picks up an empty pallet that has just been Figure 16 - Link-node arrangement, split-hold 
unloaded by the two-man team working in gang loading maritime cargo 

the port wing. The empty pallet is lifted 

out of the hold and returned to the apron. The operations of the hook are reversed during the 
next cycle. This arrangement of links is shown schematically in Figure 16. 

The results are presented as a set of time bars (Figure 17). The added amount of 
induced delay caused by the cycle-to-cycle variation is small as a result of the imbalance 
among the transporting agents and the delays caused by the work methods themselves, If the 
system were more balanced, the productivity would increase, but so too would the percentage 
of induced delay, according to Figures 6 and 7. 
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Figure 17 - Time bars, generalized loading system 
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SUMMARY 

In summary, the system of transporting material by a sequence of independently con- 
trolled movements and handlings has been idealized as a system of links and nodes, and the 
rate of flow (productivity) described as a function of six parameters. These were the number 
of links, the arrangement of links, the number of transporting agents in each link, the relative 
size of the unit of commodity handled in each link, the allowance for temporary storage at each 
node, and the mean and variance of the time required for the individual transporting agents to 
perform the basic operations of pickup, transport loaded, release, and return empty. The 
recognition that the time required to perform each operation does vary from cycle to cycle 
established the process as stochastic. Mathematical analyses of this stochastic process have 
been accomplished, providing valuable insight into the mechanism of flow of discrete entities; 
but in order to study the more complex systems which are actually used to move material, 
simulation was employed. Thus the six parameters were investigated by a series of point 
calculations, expressing the results in terms of effectiveness—the ratio of the productivity of 
the entire system to the productivity of one of the individual components. With these results it 
has been possible to describe quantitatively the phenomena which have been observed only quali- 
tatively. In addition, the calculations have revealed situations which are not intuitively obvious, 

Simulation is a powerful method of analyzing transportation networks, requiring no 
complicated mathematics. Although this paper has treated the problem of cargo handling, it is 
important to recognize that the same technique could also be used for other transportation 
systems and for other stochastic processes, such as the maintenance and service of materials- 
handling equipment. 
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MATHEMATICAL PROGRAMMING AND EVALUATION OF FREIGHT SHIPMENT SYSTEMS 


PART I - APPLICATIONS 
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This paper presents operational considerations involved in making the 
mathematical program into a working tool for management in the 
instance of crew scheduling and plant location. 











In this paper we shall present some of the results that we have achieved to date in our 
effort to apply linear programming. Our efforts along this line began about two years ago when 
staff members of the Union Railroad joined with researchers from the Graduate School of Indus- 
trial Administration at Carnegie Institute of Technology to set up a cooperative program of 
operations research. The Graduate School of Industrial Administration was anxious to use a 
complicated area, like railroading, as a testing ground for the new mathematical and scientific 
techniques of management being developed there. On our part we hoped that this testing might 
yield methods and procedures which could ultimately be installed on the railroad and which 
would raise our level of operating efficiency. After considerable exploratory discussion by the 
two staffs, the problem of scheduling crews to runs on the main line was selected as the prob- 
lem field most likely to meet the objectives of both groups. 

So that a clearer idea of what was involved in this problem may be had, it might be well 
to begin with a thumb-nail sketch of the Union Railroad facilities. The Union Railroad is a part 
of the original lake to the mills railroad which was begun in 1896 from Lake Erie to Pittsburgh, 
Pennsylvania. The Union is a Class I switching railroad, engaging in general railroad trans- 
portation and providing railroad service to many large and small plants in the Monongahela 
Valley. Our customers include such firms as the Westinghouse Electric Corporation, Linde 
Air Products Company, various mills of the United States Steel Corporation, Continental Can 
Company, and Fisher Body Division of General Motors Corporation, to name a few. We operate 
inan area in Allegheny County of approximately ten miles radius and have 500 miles of main 
line and yard track. We interchange with seven different trunk line railroads at 24 interchange 
points and employ almost 300 train crews daily. The statistics of the Interstate Commerce 
Commission show that, on a "tonnage handled basis,"' the Union ranks among the top 15 railroads. 

Our transportation operation is divided into two broad categories, namely, main line 
operations and switching operations, a classification which reflects both the nature of the work 
done by the crews and also our labor agreements. The main line work, which involves 30 to 40 
crews per day, consists of picking up cars at interchange yards with trunk line railroads and 
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delivering these cars either to other trunk lines connecting at other points or to industrial con- 
signees at other places on the road. For some customers, delivery merely means setting off a 
cut of cars at a plant or public siding. In a large majority of the plants, however, whole train 
loads of cars are placed into inbound storage yards and the customer will, at some later time, 
notify the Union Railroad when particular cars are required for unloading. Main line move- 
ment resumes when the empty cars or outbound loads are picked up in the customer's outbound 
yard and again assembled into trains for delivery to trunk lines or to some other customer on 
the railroad. Switching, as opposed to main line work, covers all car movement within a plant. 

In order to schedule main line freight trains, it is necessary to know first what cars or 
commodities are at the various interchange or plant areas and where these cars or commodi- 
ties are bound. Having this information, the Train Movement Director's office issues the 
various orders required to effect the deliveries or moves. When the staff at the Train Move- 
ment Director's office reports for work at the start of a turn, they will find a list of the loco- 
motives of various sizes and capabilities available at approximately 21 different starting points. 
In addition, they will also know how many crews have been hired for that particular turn to 
operate the locomotives and where these crews will report for work. With this information, 
specific assignments of work to be performed by each crew during part or ail of the turn are 
drawn up. The Train Movement Director has a wide variety of possible assignments to choose 
from and also a large number of restrictions on his choice. Frequently he must plan for 
pusher service, as when a crew must meet another crew along the route to assist it up a par- 
ticularly steep grade. He must decide whether a crew should gather together a solid train load 
of cars all bound for one destination and haul them through to the destination; or, in the belief 
that lower classification costs and longer trains may offset losses in cross-hauling, he may 
decide to have the crew pick up mixed lots of cars at origin stations and carry them to an inter- 
mediate classification yard to be sorted and regrouped for the next stage of their journey to 
their ultimate destination. Crews must begin and end their shift at the same place, and since 
dead-heading is normally not allowable, this restricts the work assignment. Crew performance 
varies with a given shift, and since crews are hired on a strict seniority basis, assignments 
are planned with the capabilities of the crews in mind. 

This then is our main line scheduling operation and some of the factors affecting sched- 
uling decisions. How well are we doing this job? How could we do it better? Because of the 
complexity of the operation, these were not easy questions to answer by conventional methods. 
That is why we thought that an operations research approach or mathematical simulation of the 
problem might be useful. 

In experimenting with various types of mathematical formulations, the model develop- 
ment followed a pattern which would simulate our own scheduling operation. We recognized 
from the outset that it would be impossible to try to deal with all phases of the operation at 
once. It would be necessary first to make many simplifications of the scheduling problem and 
then to test the method developed step by step, bringing in as we went along the additional 
features necessary to enable us to utilize the method in practical day-to-day operation. A 
number of different formulations were tried, and when Doctors Charnes and Miller had a for- 
mulation which seemed to capture the essential parts of our scheduling problem, we then began 
a testing program to determine the feasibility of the formulation, as well as to determine the 
refinements which had to be made to it. We devised a tentative coding system which consolidated 
the stations on the road into seven basic points. (Refer to Figure 1.) Eventually we shall prob- 
ably have to expand the number of points, but this division was fine enough for testing the 
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Figure 1 - Coding system of physical area serviced by 
union railroad 


approach and yet lead to a problem still small enough so that we could conveniently 
compute the program by hand. 

An idea of what was involved in the computation is given in Figure 2, which shows a 
typical working tableau. On this chart we list in the left-hand column the routes to which 
crews and engines can be assigned. Actually, there are over one hundred combinations, and 
from these we selected approximately eighty feasible routes to use. Since our crews are hired 
on an eight-hour-shift basis, we have assumed that no more than three places or points could 
enter in to any one assignment for a crew. In practice, we find that today we get very few crew 
assignments where more than three different points are reached by a crew in any given shift. 
These route assignment possibilities will be expanded such that a crew can traverse more than 
atwo- or three-leg assignment in an eight-hour shift if methods are improved or if we get 
better performance from the crews. 

In the bottom horizontal row, we list the requirements in terms of train loads of mate- 
rial that must be moved between various points on the railroad. This information is available 
tous at the beginning of the shift. 

In the unit-cost column, we indicate the cost for a crew unit of service on the route 
involved. For example, in the cost column, we list 1 for route 1,2. This could indicate the 


relative cost for providing crew service or could indicate the total crew cost for an eight-hour 
turn. 
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It should be noted that, for example, in moving a train load of material between points 1 
and 2, one merely has to look in column 1-2 to see that service can be provided by utilizing any 
one of five different routes. Likewise, the same would be true for providing service between 
any other points. Other legs of routes can be used to provide the service. 

The problem now becomes one of deciding which combination of routes should be used 
to provide low-cost movement. If a systematic procedure is followed in solving this over-all 
problem, a complete scheduling program will emerge which will indicate the route assignment 
for each crew; and the integrated assignments will insure that this is at minimum cost, which 
of course implies that we have the minimum number of crews making these movements. 

For example, we have five train loads of material to move from point 1 to 2. Our 
mathematical solution tells us to schedule the movement of the trains as follows: 

1. Have two train loads moved by scheduling 2 crews on route 1,2. Of course, this is 
only part of the work for this crew. It will also move trains from 2 to 1 as part of this 
assignment. 

2. Another crew, assigned to route 1, 2,6 will pick up another train load of material at 
point 1 and move it to 2. This also is only part of the assignment, for it too will pick up a train 
at point 2 and move it to point 6 and a train at 6 and move it to point 1. 

3. To complete the shipments between points 1 and 2, our plan shows that two crews 
should be assigned to route 1,2,7. In this case two train loads of material will be moved 
between points 1 and 2 as part of this assignment and, in addition, this crew will carry loads 
between 2-7 and 7-1. It should be noted that there is only one train load of material to be moved 
from 7 to 1, thus one of the crews will go light between 7 and 1. Significantly, we are schedul- 
ing light moves where it costs us the least. 

4. It can be seen that all the shipments between the same points are made by assign- 
ment to different crews. Where there is nothing to move, the crew will travel light, and this is 
clearly indicated in the light or surplus row. 

Not only does this computational method tell us what the best combination of crews is 
and when we have found it, but it generates an additional set of numbers called evaluators. 
These evaluators, shown in the horizontal row at the top of the chart, have potentially great 
significance for our management. For one thing, they tell us how many more crews we would 
have to hire if we had one more train load of traffic on the leg in question. 

Another use: Suppose we were considering some location on the railroad for a plant 
location; this is especially helpful when we are encouraging some company to locate in an area 
serviced by our railroad. By computing the evaluators, we can determine where such a loca- 
tion could be serviced at least cost to the railroad. An evaluator of zero would indicate that 
additional service could be provided to this area without the railroad incurring any additional 
crew cost. Thus we might suggest to the potential customer that favorable rates would apply if 
this location were selected by him. 

Today, railroad management is especially inte.ested in seeing some Congressional 
action on the 'Weeks' Report." One of the important points in this report is the "setting" of 
Compensatory tariffs on a competitive "cost basis" rather than on an industry-wide negotiated 
basis. The relative cost to make freight shipments between various places is one of the 
"dividends" of this mathematical solution and should be useful in strengthening our position to 
secure pricing relief in accordance with the recommendations contained in the Weeks' Report. 

The testing consisted of verifying the mathematical formulation by selecting several 
days' actual main line movements which had been made. The actual number of crews employed 
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to move the train loads of material, and also the light moves, were recorded. The minimum 
number of crews necessary to move the same number of train loads of material that were 
actually moved was then determined by going through the computational routine in the mathe- 
matical formulation. We found that the optimal routing and freight scheduling was feasible in 
all cases. A comparison of the results obtained by the new method with those of actual operat- 
ing experience showed that there are also potentially large savings to be made, possibly 10 to 
20 percent of the crew costs if the mathematical approach to crew scheduling is used. This 
large margin for cost reduction should not be considered, however, as reflecting on the com- 
petence of our scheduling people. The trouble is that there are literally thousands of ways of 
scheduling crews for any given shift, even after allowing for all the many restrictions govern- 
ing the crews, traffic, and facilities. Human beings relying solely on judgment, no matter how 
expert, can only examine a few of the more likely-looking prospects in the limited time avail- 
able. The mathematical approach, by contrast, makes it possible to examine in effect all the 
combinations of runs which are feasible and to select from this mass of possibilities the one 
or more which are best as measured by some tangible standard, such as least cost to the 
railroad. 

In applying the mathematical routine, we used regular clerks, so we’can conclude that 
the new method was found to be technically feasible in the sense that it could conceivably be 
used by the kind of personnel we are likely to have. With it they could succeed in translating 
shipping requirements into our type of schedule of crew runs between points. The method 
would definitely not require the services of trained mathematicians but neither would it be com- 
pletely a mechanical clerical job. A considerable amount of judgment must be applied to get it 
from the numerical results to the concrete program, but we believe that the judgment that is 
needed is the kind that our scheduling people should easily be able to make. 

As for the future, we hope to continue toward our goal of making the mathematical pro- 
gram into a working tool of management. Before this can be done, it will probably be necessary 
for us to add additional restrictions of various kinds to the basic model with which we have been 
working. Just which restrictions we will have to add directly and which we can leave to judg- 
ment can be settled only after further computations and study of the patterns which emerge. 
For this stage of the work, where so much computation will be required, we shall probably have 
to have access to a high-speed electronic computer. We don't mean to imply that it will be 
absolutely necessary that we have a high-speed electronic computer to effectively utilize this 
mathematical formulation in our future application work, for we hope to develop some other 
routine by which these computations can be done even more quickly. 

It may also prove to be the case that after we have studied enough programs we will 
begin to see more clearly whether there are any systematic sources of inefficiency in our 
present scheduling arrangements. These, of course, can be corrected directly, even in advance 


of being able to turn over to the Train Movement Director a perfected model to do the whole 
scheduling job. 


In conclusion, we at the Union Railroad are very optimistic about the progress of our 
work and the mathematical formulation of freight shipment scheduling. We believe that what 
we have here is a basic development which probably will find wide application in railroad indus- 
try and undoubtedly in the transportation industry as a whole. 
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Models are presented which permit simultaneous analysis of crew 
assignment, routing, and switching and their interaction, particularly 
in consideration of circuitous routing versus expansion of classification 
and switching facilities. 











I. INTRODUCTION 

The previous paper by Mr. Soyster was concerned with the "application" side of pro- 
gramming freight movements on the Union Railroad. It described the specific setting of the 
operation; the progress to date in developing and testing a mathematical model of the problem; 
some of the ways in which that particular model can be used at both the operating and planning 
levels; and something of the plans for pushing ahead with the work and bringing it to full fruition. 
The concern in this paper is with certain structural aspects of the model* and with an example 
of how the model can be adapted and extended to cover a wider class of transportation problems. 


I. THE JOINT-COST PROBLEM 

The so-called "transportation" or "distribution" model has been applied successfully to 
some facets of the railway problem, notably the routing of empty boxcars. In view of the great 
computational simplicity of that model, it was natural to begin by attempting to extend the dis- 
tribution model to car movements between origin and destination points on the Union Railroad. 
It quickly became apparent that the distribution model, as it stood, could not be so extended. 
There were several reasons why this proved to be the case, but the most serious was that the 
cost of moving a unit of traffic—even a gross unit like a full trainload—between any two points 
was unknown. By unknown is meant not simply that accurate cost information was hard to get 
from the accounting records, although that was certainly true. The difficulty was rather that 
each train crew, having completed a delivery between two points, could then go on to make one 
or two or more deliveries between points during the eight-hour shift and possibly to do some 
other productive work as well, such as switching in one of the yards. In fact, each train crew 
had to be so routed as to return to its starting point within the shift. Hence a delivery between 
apair of points was never more than half, and normally would be a good deal less than half, the 





*For a discussion of the concept of ''modelstructure'' and its significance, see Reference [1]. 
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total potential work-load of a crew. Under these conditions the "cost" of any particular ship- 
ment between points cannot be determined without knowing the best possible routing pattern 
for all the shipments. 

This problem is essentially an extreme case of what economists refer to as joint pro- 
duction and, hence, joint costs. In the literature of economics, railways have always been the 
classic textbook examples of joint production, and if one looks at actual railroad operations— 
both at the Union Railroad and elsewhere—one can easily see many instances of the phenom- 
enon. Something like the crew-return-to-base requirement is found on every railroad; and the 
same general sorts of looping restrictions apply to rolling stock as well. Obvious counter- 
parts to these requirements exist in other forms of transportation. 

To capture the essentials of this problem, a model was needed which could do two 
things: First, at the operating level, it had to provide a method for determining optimal work 
patterns when the basic productive unit could be assigned to several distinct productive tasks 
during a single scheduling period. In the specific context of the Union Railroad this meant a 
model which could determine the best routing assignments for train crews making deliveries 
between pairs of points on the network, and which could make the assigaments subject, among 
other things, to the constraint that each crew be routed back to its starting point within the 
eight-hour turn. Second, the model had to provide a method for costing or evaluating the vari- 
ous components of the individual work assignments. Such information, though it might not be 
necessary for the purely operating part of the process, is indispensable for facilities planning 
and for related types of management planning in situations of joint production. To turn again to 
the context of the Union Railroad, this meant a model which could tell management how much it 
really cost to make another shipment between any pair of stations, in the sense of allowing for 
ail indirect repercussions of the added traffic throughout the system as a whole. 

A simple and revealing model meeting these tests is the following. To make the very 
striking mathematical structure of this model stand out sharply and to facilitate comparison 
and contrast with other model types, it will be helpful to proceed directly to a simplex tableau 
of the problem (Table 1). 

In the requirements vector, Py are entered the minimal shipment requirements in 
units of trainloads between pairs of origin and destination points during the scheduled period. 
On the Union Railroad the actual amounts to be entered here on any shift would be determined 
in the last analysis by customer calls for cars, scheduled arrivals or departures of road trains 
at interchange yards, ore dumper schedulers, and the various other factors affecting car move- 
ments as described by Mr. Soyster. 

The next set of vectors, labeled R(1,2), R(1,3)...R(i,j)...R(m,n); R(1,2,i)...R(i, m, n), 
represent the various routings to which crews and engines can be assigned. The route R(1,2), 
for example, is a two-leg round-trip route running from point 1 to point 2 and back to point 1. 
Each crew assigned to that route could carry one trainload outbound from 1 to 2, hence a coef- 
ficient of +1 is entered in that column in the row corresponding to shipment requirements 1 and 
2. On the return trip, a trainload could be carried from 2 to 1, so a +1 is entered in that row 
also. Thus the return-to-starting-point constraint in this formulation is incorporated directly 


into the matrix by entering the +1's in the rows in such a way that the destination on the last leg 
is the same as the origin of the first leg. 





*In the tableau only two-leg and three-leg possibilities have actually been set down, since, on the 
Union R.R., runs involving more than three legs are rarely feasible within an eight-hour shift. 
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TABLE 1 
Structural Tableau of Train-Routing Model 





Ci; + Gp» C)3 eee Gj eee Coe eee Choi eee Cows M M eee 
P, R(1,2) R(1,3)...R(i,j)...R(m,n)...R(1,2,i)...R(i,n,m) S(1-2) S(2-1)...S(n-m)...(1-2) (2-1) 





1 1 -1 1 
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Above the routes, in row labeled C; are entered the costs of assigning a single crew- 
engine package to the route in question. These costs may be stated in a variety of ways depend- 
ing on the specific purpose of the moment. For example, they may be taken as the standard 
crew and engine expense; or as the expected expense allowing for the fact that on longer runs 
there is a greater possibility of overtime. In other cases, and this was true at the Union, it is 
sufficient to work in terms of minimizing the number of crews assigned, which means treating 
the unit cost of each route as 1. 


The next set of column vectors after the routes in the columns labeled S(1-2), S(2-1)... 
§(m-n) are the slack vectors. Since this is a minimization problem and since all the inequalities 
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are in the form of requirements that shipments be equal to or greater than the minimums 
specified, these slack vectors represent overfulfillment of the requirements. Hence, they 
require an entry of -1 in the appropriate row. In terms of the railroad operation they repre- 
sent "light moves," i.e., movements of crews, and engines without cars. 


The last set of vectors, labeled (1-2), (2-1)...(n-m), are artificial unit vectors neces- 
sary to form the initial basis in a standard simplex tableau, They may be referred to as "leg 
vectors,"' because they represent, as it were, the dismembered legs out of which feasible 
routes—which must contain at least two legs—can be built up. That is, any feasible route can 
always be expressed as a simple vector sum of corresponding leg vectors, a fact which plays a 
key role in the special computation routines developed for obtaining optimal programs. 

In analytical terms, let 
@ Ryw.in....b. .) the i*® routing loop whose legs Ly hy yeh 

j,k" k,1 P,4 connect the points j,k...q, and which may be 
o='2.. traversed from j through q and back to j within 
some specified time period 


P,q 


X; fie the number of loaded trains sent from j to k over 
j,k the th routing loop R; AL L. = 
g,h°°°"j,k°** "p,q 


the number of light trains sent from j to k over 


g, n°? j, °° BL e 


the total number of trains, loaded and light, assigned 


sth 
to any leg of the i— route R, 
(Lg nee e Ly ieee Ly q) 


the total of all trains, loaded and light, shipped from 
j to k over all routes containing a leg from j to k 


the cost of assigning one train and one crew to the 
jth route 


the total number of loaded trains required to be 
shipped from point j to point k during the specified 
time period. 


We know from the requirement that all trains and crews must finish their runs at the 


same point at which they started that for any route R. 
i/(L, g°+++ Lng) 


X, + 3 
i/Ly 9 


X, 
/L, r+1 
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Hence, we may express (8) more compactly as: 


(9) X, for all Lei and for all i. 


i/ 
L r,r+l 
With these definitions it is possible to formulate the routing problem in either of two 
equivalent forms: 
Formulation I (Leg Formulation) 





(10) 


Subject to: 
m 
(11) Le Xi y 2 bj, for all pairs of points j, k; (j * ). 


(12) X, > X, for all Lei and for all i. 
. Ly 4 


(13) X,, X 0. 


: > 
V/Lj_ 


Formulation IT (Route Formulation) 





(14) 


Subject to: 


(15) »: X; 2 Dik for all pairs of points j, k; (j ~ k). 


ie(hy, 


(16) X,2 0. 


The second or route formulation is the one which was employed in constructing the 
simplex tableau of Table 1. 

That these are equivalent formulations can be seen by supposing that we have an optimal 
solution to I. The X; associated with it must also necessarily satisfy II, and hence the optimal 
value of I must be equal to or greater than the optimal value of II. Suppose further that we 


have an optimal solution to I, say X* . Then we can define some x4 /L = x} , and all other 
r,r+l1l 
sin = 0. We would, therefore, have a solution to I with the same value of the function to 
J 
be minimized. 
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The structure of the coefficients in this model has a strong family resemblance to the 
familiar distribution model. There are important differences, however, which are sufficient to 
rule out the possibility of applying one of the many standard computational routines developed 
for the distribution model. But while the distribution model routines will not serve, the fact 
that the structure of the routing model is still a fairly simple one, as structures go, does per- 
mit the employment of techniques which at least are a good deal more efficient than a straight- 
forward simplex computation. A detailed description of the routine which was developed to 
exploit the +1, -1, 0 structures of the routing model is given in Reference [3]. The general 
outlines and underlying logic of the routine, however, can readily be sketched in. 

The routine is based on the so-called modified simplex method of Charnes and Lemke (2), 
That method relies on the fact that the inverse of the basis at any stage of a simplex calculation 
provides all the information necessary for making alterations in the basis to improve the pro- 
gram. That is, if we take the following augmented inverse: 


Bp! o 


w 1 


where w is the m x 1 row vector representing the evaluators or indirect costs of the column 
vectors composing the inverse, that is, w= c we , and multiply it by the original value of 
— 


any vector not in the basis “ , we obtain the vector “i , where the zi- cj represents 
~C; Z5 - Cj 

the net change in the functional from bringing the vector into the basis; and x are the coeffi- 

cients expressing the vector in terms of the vectors currently in the basis. The rest of the 

calculation proceeds according to standard simplex procedure. 

In a simplex tableau, the inverse of the basis at any stage appears as the entries under 
the unit vectors which formed the initial basis. In the routing model, these unit vectors are the 
leg vectors, the fragments out of which feasible routings are built up by summation subject to 
the looping constraint. Hence the tedious process of matrix multiplication in the ordinary 
modified simplex routine can be replaced by simple vector addition. Any vector not in the 
basis can be expressed in terms of the current basis vectors merely by summing the appro- 
priate columns of the inverse. Furthermore, since the +1, -1, 0 structure of the original 
matrix normally carries over to subsequent tableaus, the process of adjusting both the program 
and the basis inverse can also be reduced to a matter of vector addition and subtraction rather 
than the laborious divisions and multiplications of the normal simplex routine. * 

Not only can this formulation of the joint routing problem be exploited to provide effi- 
cient routines for generating optimal work assignments, but it also provides an answer to the 
second part of the joint production problem, namely, the costing of the various components of 
the workload. These costs or evaluations are given directly by the optimal values of the vari- 
ables in the dual of the routing problem. That is, each dual variable represents the marginal 
cost of moving another trainload between the points in questions, mutatis mutandis. There is, 
of course, no need to solve a new problem to get these marginal costs; they emerge automati- 
cally as by-products in the solution of the direct assignment problem itself. 








*Further time-saving comes from the ability to start with a basis fairly well along, since its 
inverse can be generated relatively easily. 
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[J]. PROGRAMMING AND EVALUATION: AN ILLUSTRATION 

The basic model described in the previous section is intended to capture the essence of 
the crew assignment problem. To apply the model on a day-to-day basis, or to use it at the 
planning level, it is necessary to include additional features of the operation, particularly those 
relating to timing and to the fact that future shipping requirements are uncertain. Extensions 
along these lines are currently under development, and an account of them will appear in due 
course in the future. Rather than dwell on these, it would seem to be more in keeping with the 
announced objectives of this conference to explore possible extensions of the routing model to 
programming and evaluation in other transportation contexts. 

One such extension, which may be of interest to the trunkline people present, concerns 

certain aspects of the switching problem. Consider, for example, the simple three-point rail- 
way system of Figure 1, where each of the 
three points represents a major complex aie ee 
of industrial and classificationyards. Cars 1,2*.5 
originating at the industrial yards within, hg . 
say, point 1 and destined for industrial or ae 
interchange yards in, say, point 3 must 
somewhere along the line be classified 
according to the particular industrial yard 
at point 3 for which the cars are destined. * 
This sort of classification is made presum- 
ally so as to simplify the actual delivery 
operation at the receiving end. Under the 
conditions of the problem it is assumed 
that the actual classification may be accom- 
plished in any of three alternative ways: 





(1) The classification may be made 
atthe origin point 1, that is, by "blocking" 
the train in station order in the outbound yard and routing that blocked train directly to 3. 

(2) Or, the various cars can be sent out from 1 to 3 in a mixed train and the cars 
switched into suitable delivery units by crews of the classification yard at 3. 

(3) Or, finally, a mixed train of three-bound cars can be routed from point 1 directly 
via the hump yard, which is assumed to be located at point 2. At point 2 these cars will be 
"blocked" in station order and sent on to their ultimate destinations at point 3. 

The actual amounts which must be shipped between the various points during the sched- 
we period are indicated on the map by the numbers with directional arrows. (Thus, there are 
tine trainloads of cars to be shipped from point 1 to point 3; five trains from 3 back to 1, and 
s00n.) The average switching cost (S,) and the maximum switching capacity (K,) of each yard 
is shown next to each point. These capacities may be thought of either as total physical proc- 
essing capacity or as the capacity available after other yard work is scheduled. The yard at 
point 2 is assumed to be a hump yard, with the very large capacity of 24 trains a period and the 
lowest switching costs of .1. Costs are higher and capacities lower at the other two yards.** 


—_—_—_ 

*Other ''processes,'' such as inspection, can readily be substituted for the switching of the 
example without chauging the formal outlines of the problem. 

"Needless to say, no particular significance should be attached to the absolute values of these 
humbers. They were chosen pretty much arbitrarily, though with some effort, to get relative 
Costs into a reasonably meaningful pattern. 
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In addition to switching costs, there are, of course, train movement costs, and we shall 
assume these are of two kinds: 

(1) First, running costs, consisting of fuel consumed plus the wear and tear on the 
engines. The actual values used for running costs between points. are shown in parentheses on 
each leg. For example, on the leg from 1 to 2 the running cost is .2 (this being an uphill run), 
whereas on the down side from 2 to 1 the cost is less, .1. An allowance might also be made 
for the difference in cost of an engine moving loaded compared with one moving light, but for 
simplicity this difference was assumed to be negligible. 

(2) The second movement cost is, of course, the train crew cost. Road crews on this 
railroad are assumed to work under agreements similar to those on the Union Railroad—that is, 
they must be routed back to their starting point within the scheduling period. This means, as 
was discussed earlier, that crew cost cannot be allocated to the particular legs of a run. All 
we know is the cost of assigning a crew to one of the routing loops. For this problem, to keep 
matters simple, we will consider only two-leg routes, and we will use the same figure 0.5 for 
crew costs on each of the loops.* 

Given this pattern of requirements, capacities, and costs, we want to know such things 
as the following: 

(1) At the operating level we want to know how to route and where to classify each train- 
load in such a way as to minimize over-all switching, running to labor cost. We especially want 
to find out to what extent it will pay to route trains circuitously via the low-cost hump yard, that 
is, to substitute road expense for switching expense. 

(2) At the planning level, we want to know which of the yards is the real bottleneck in 
the system. That is, if we have a limited amount of funds to invest in expanding yard capacities, 
where will an extra unit of switching capacity give us the greatest return for our money? Per- 
haps, under current conditions, the more interesting question might be not so much to find the 
bottleneck as to find the empty bottle. That is, assuming we want to reduce the number of yards 
we operate, how much does it cost in terms of extra running, crew, and switching expense for 
us to shut down any particular yard? And what will be the implications of closing down this 
yard on traffic flows and switching levels throughout the rest of the system? 

(3) Also at the planning level, we would want to know something about the cost of ship- 
ment between the various points. Once we have some idea of which are the high-cost and which 
are the low cost legs, we have some idea of where the traffic department should be concentrat- 
ing its sales promotion so as to have the greatest chance of bringing in new traffic which can 
be stuffed into the slack in the system. 

It can readily be seen that this problem, though simplified in many respects, is now a 
good deal more complicated than the simple joint-routing problem. Yet it turns out that, even 
with this additional detail, the mathematical structure of the model remains quite simple, and 
the structure of the model preserves those features which were exploited in the basic model to 
provide a rapid computation routine.** 





*There is nothing crucial about the up-and-back-within-the-shift requirement. It is entirely 
feasible to allow for such things as layovers or crew transfer points along the legs. All that 
is important is that the crew restrictions involve joint production, that is, that crew cost 
cannot be allocated to particular legs. 

**For computational purposes the problem is more easily handled in a "leg formulation" (see 
p. above), with additional constraints covering switching capacity and insuring that inputs 
to the hump yard on indirect routings equal outputs from the yard. Considerable progress has 
recently been made in developing special routines for this model and, in fact, for a much 
wider class of related problems. An account of these will be forthcoming in the near future. 
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ANALYSIS 


The numerical solution of the (direct) problem is as follows: 

(1) Assign six crews to the routing loop 2, 3 and switch at 2 all traffic both from 2 to 3 
and 3 to 2 (making nine trains to be switched at 2). 

(2) Assign four crews to the routing loop 1, 2 and again switch all traffic both ways at 2 
(making an additional five trains to be switched at 2). 

(3) Assign six crews to the routing loop 1,3. Since there are nine trainloads from 1 to 
3 and only six crews on the route, the three excess trainloads are to be routed "around the 
horn" and switched at 2. This makes total switching at 2, seventeen trainloads. There are 
eleven trainloads left on the 1, 3 loop to be switched at either 3 or 1. The model tells us to 
switch nine of these at 1 and two of these at 3, but beyond this it does not make any difference 
which particular shipments are sent to the one yard and which to the other. Since we have six 
crews on the loop 1, 3 and only five trains moving from 3 to 1, we have one light train on this leg. 

With respect to routing, then, it was profitable to engage in some circuitous routing 
"around the horn."" But as the numbers stand, this was profitable not so much because switch- 
ing costs at the hump yard are only 1/3 those in the other two yards, but rather because traffic 
flow into and out of point 2 is unbalanced. That is, on direct up and back routing we would have 
three crews traveling light from 1 to 2 and another three traveling light from 2 to 3. Since we 
have these crews available, routing three trainloads of 1-3 traffic via 2 enables us to eliminate 
three crews from the direct 1-3 run. Notice that it must be 1-3 traffic that gets the circuitous 
routing not 3-1 traffic; if any traffic is sent the other way, the extra crew cost will much more 
than offset the saving in switching costs. 

As for the cost of shipment over each of the legs, as given by the solution to the dual 
problem, there is a considerable range of variation.* On some legs, the solution tells us that 
an extra train could be handled for just the switching and running costs, no more crews being 
required (e.g., a trainload from 3 to 1). At the other extreme, an additional trainload will, in 
some cases, mean that the expense of an additional train crew will have to be incurred in addi- 
tion to the running and switching costs (1 to 3, for example). In between are some cases where 
an extra trainload will require some extra crew expense, but less than proportionately so 
because of certain offsetting adjustments which can be made elsewhere in the system. 

Insofar as 1-2, 2-1 traffic and 1-3, 3-1 traffic is concerned, the relative costs on each 
leg, at least, are in accord with orthodox expectations, namely, the cost is higher in the heavy 
or main directions (1-3 and 2-1) than in the light or back-haul directions (3-1 and 1-2). With 
respect to 2-3, 3-2 traffic, however, the situation is strikingly different. Here, taking on an- 
other trainload in what looks like light direction 2-3 will cost the railroad considerably more 
than adding another train in what looks like the heavy direction 3-2. This peculiar result is 
explained by the fact that adding one more train load from 3 to 2, which means another poten- 
tially light run back from 2 to 3, permits routing one more trainload of 1-3 traffic around the 
horn. One more crew is added on the 1-2 loop, thereby, but one can be dropped from the 1-3 loop. 
The net result of all these changes is that part of the extra crew cost is offset by lower switch- 
ing costs. 





*The costs between pairs of points were as follows: 
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In the matter of bottleneck yards and potential areas for expansion, it turns out that 
none of the yards merit an increase in capacity at these traffic levels. This is obvious enough 
in the case of the yards at 2 and 3, which clearly have excess capacity as is. But it also is 
true of the yard at 1, which is apparently busily operating away at full capacity. If by some 
chance more traffic should arise between 1 and 3, the extra switching load can be handled at 3 
at no greater cost than at 1. What counts in this case is the combined capacity at 1 and 3; and 
as long as either one has excess capacity, neither is worth expanding. 

While there is no basis for choice between yards 1 and 3 when it comes to expansion, it 
makes a good deal of difference which one is closed down. Closing down the yard at 3, for 
example, will increase aggregate running, switching, and crew costs in this system by less 
than 3 percent—an amount small enough so as likely to be offset by the capital recovery from 
liquidating the facilities. Closing down yard 1, however, will increase costs by about 15 per- 
cent. The reason for this difference lies not in the yards themselves but in the pattern of traf- 
fic and the adjustments to schedules which must be made to compensate for the decrease in 
available capacity—and this even though the system as a whole has, and will still have, excess 
capacity. When the yard at 3 is closed down, the pattern of traffic is such that its workload can 
be absorbed by compensating movements around the horn—one train from 1 to 3 and one train 
from 3 to 1. Thus the number of extra road crews required is kept to the minimum. In the 
case of yard 1, however, no such fortunate result occurs; redistributing its work leads toa 
substantial imbalance in the flows around the system. 

Such, then, is one possible type of extension of the basic joint-routing model. Obviously, 
since the numbers are artificial and since many phases of actual routing and switching opera- 
tions were deliberately ignored, we cannot draw any action-conclusions from the numerical 
results themselves. The model, however, simple as it is, is not totally academic. The model 
illustrates, among other things, the dangers in the currently fashionable view that the main line 
routing operation is in fine shape and that the problem is in the yards. This is not to say that 
the mechanical aspects of the switching operation are not well worth studying or that potentially 
large savings cannot be made from consolidating switching facilities. It is easy to forget, how- 
ever, when becoming immersed in these individual yard studies, that a switching facility is not 
an isolated property. It is a part of a much larger and finely interlocking system. For effec- 
tive top-level planning of yard-location policy under these conditions, it is essential to have an 
approach capable of dealing with the system as a whole, routing as well as switching. 





*That is, the dual variables in all three cases were zero. 
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NUMERICAL SOLUTION OF THE GENERALIZED TRANSPORTATION PROBLEM 


William Prager 
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A numerical "saturation'' technique is presented for solution of a 
generalization of Hitchcock's distribution problem, 











1. INTRODUCTION 

A generalization of Hitchcock's [1] transportation problem was discussed in a recent 
paper [2]. This generalization is concerned with the steady flow of a homogeneous product 
from several production centers (plants) to numerous consumption centers (markets). A plant 
and a market specify a shipping route, and the cost of maintaining the unit rate of flow through 
the unit time is known for each shipping route. This specific shipping cost is supposed to be 
independent of the flow rate along the considered route. Furthermore, a route may have a given 
capacity, i.e., a maximum flow rate that it can accommodate. 

With each center, there is associated the rate of production or consumption and the 
price that has to be paid for producing or consuming the product at a unit rate. This specific 
price is supposed to be independent of the rate of production or consumption. The prices at a 
market and a plant establish the price differential for the shipping route connecting this market 
to this plant. 

A steady flow is stipulated to satisfy the following conditions: 

(i) No route carries a negative flow (i.e., a flow from market to plant) or a flow in 
excess of its capacity, and the total amount shipped from each production center per unit of 
time or received at each consumption center per unit of time equals the rate of production or 
consumption rate at this center. 

(ii) The price differential for a given shipping route is not larger than, equal to, or not 
smaller than the specific shipping cost for this route, according to whether the route carries 
zero flow, positive flow below capacity, or capacity flow. 

With this terminology, the generalized transportation problem can be stated as follows: 
given either the price or the rate of production or consumption at each center, to determine 
the other prices and rates at the centers and the flow rates along the routes in such a manner 
that the conditions (i) and (ii) are fulfilled. 

The uniqueness of the solution of this problem was investigated in the paper already 
mentioned [2], and extremum properties of the solutions were established. The present paper 
describes a method for the numerical solution of problems of this kind. Ford and Fulkerson [3] 
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have given a similar method for the numerical solution of the Hitchcock problem. The method 


presented in this paper is a natural extension of their method to the generalized transportation 
problem. 


2. TYPICAL PROBLEM 

An example of a generalized transportation problem is provided by the following situa- 
tion: the rate of production or consumption is given at each plant or market; the total rate of 
production of which the plants are capable does not equal the rate at which the markets can 
absorb the product, and a financial penalty is assessed at each market for the unit rate of un- 
satisfied demand. 

To bring this problem into the class described in the preceding section, it is necessary 
to introduce an additional production center, called the store, and an additional consumption 
center, called the dump. The store fills the otherwise unsatisfied demands, and the dump 
absorbs the production excess resulting from the possibility that it may be cheaper to pay the 
penalty for not satisfying some demands than to satisfy them at all cost. Since amounts 
"shipped to the dump" are in reality not produced at all, the specific shipping cost from any 
plant to the dump is zero, and the price at the dump is zero. Moreover, the price at the store 
can be set equal to zero, provided that the specific cost of shipping from the store to a market 
is taken equal to the penalty for the unit rate of unsatisfied demand at this market. Finally, 
the part of the infinite capacity of the store that is not used to cover unsatisfied demands must 
be dumped. Accordingly, the specific cost of shipping from the store to the dump must be 
taken as zero. 

In the earlier paper [1], a system of flow rates, production rates, and consumption 
rates that satisfy the conditions (i) above was called kinematically admissible, and the adjusted 
cost of shipping was defined as the difference between the total cost of all shipping operations 
and the income from the operations at those centers where the price is prescribed. Moreover, 
it was shown that all solutions of a generalized transportation problem have the same adjusted 
cost of shipping, which is smaller than the adjusted cost of shipping associated with any kine- 
matically admissible system of rates of flow, production, and consumption that do not correspond 
to a solution. 

For the problem considered here, the prices are prescribed only at the store and the 
dump, and these prices are zero. The adjusted cost of shipping therefore reduces to the total 
cost of shipping from the plants and the store to the markets, i.e., the actual cost of shipping 
from the plants to the markets augmented by the penalties for leaving demands unsatisfied. 
Any solution to the considered problem therefore minimizes this sum. Our interest in the 
solutions stems from this minimum property. 

In the following three sections, the essential steps of the proposed method of solving 
this problem are illustrated by means of a numerical example, which involves three plants and 
four markets. The heavy lines in each of the Tables 1 through 5 bound a rectangle of twenty 
cells. These are arranged in four rows (corresponding to the plants A, B, and C and the store 
S) and five columns (corresponding to the markets a, b, c, and d and the dump m). The num- 
bers in parentheses to the right of the plant labels or below the market labels indicate the given 
maximum rates of production or consumption of which these plant and markets are capable. 
The maximum total rate of consumption (31) exceeds the maximum total rate of production (29) 
by 2, so that the rate of unsatisfied demand must be at least 2. It will be seen, however, that 
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the optimum program of production and shipping (Table 4) does not use the full productive 
capacity of plant A and involves a rate of unsatisfied demand totalling 4, the rate of unsatisfied 
demand at each of the markets a and d being 2. 

Each cell in Tables 1 through 5 corresponds to a shipping route. The numbers in the 
upper left, upper right, and lower left corners of a cell, respectively, represent the given 
specific shipping cost, the given capacity, and the flow rate for the corresponding route. When 
the upper right corner of a cell is empty, the capacity of the corresponding route is supposed 
to be so ample that it does not impose any restrictions on the shipping program. When the 
lower right corner of a cell is empty, the shipping program represented by the table does not 
involve any flow along the corresponding route. 

The numbers immediately to the left of the left-hand heavy line and immediately above 
the upper heavy line are the prices at the various centers of production and consumption. In 
particular, the prices at the store and at the dump are zero in each of the Tables 1 through 5. 

The procedure described in the following consists of a single application of Step 0 
("getting started") followed by repeated applications of Step 1 ("saturating given routes") and 
Step 2 ("raising prices"). At a generic phase of the procedure, a program of production and 
shipping is obtained that constitutes a solution to a modified problem for which the rates of 
consumption at some markets are lower than for the given problem. Step 1 raises these rates 
of consumption by using routes already available under the price system of the modified prob- 
lem. When no further progress can be achieved in this manner, Step 2 is applied; it brings 
new routes into play by raising prices in a manner that does not eliminate routes already used. 
The procedure ends when the consumption rates of the given problem have been attained. 


3, STEP 1: SATURATING GIVEN ROUTES 
Table 1 shows the solution to the modified problem obtained at a typical phase. Any 
cell for which the price differential equals or exceeds the specific shipping cost is marked by 
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a circle. The conditions (ii) in Section 1 require that a cell carries capacity flow if the price 
differential exceeds the specific shipping cost; the circles in cells of this type are marked by 
crosses. A cell may carry capacity flow even though its price differential only equals the 
specific shipping cost; the circles of cells of this type are marked by a slanting line. 

While the total flows in columns b, d, and m of Table 1 match the given consumption 
rates at these markets, the total flows in columns a and c fall short of the given consumption 
rates of these markets by the amounts noted in an additional row M below the lower heavy line, 
Table 1 thus indicates the solution to the modified problem for which the consumption rate at 
A is 9 and that at C is 5. 

Inspection of Table 1 readily reveals a way in which the consumption rate at C can be 
raised while the prices are kept constant; alternatingly decrease and increase the flow rates 
in the cells of the path Am, Ab, Cb, Cc, and Mc by the same amount q. Except for row M and 
column m, this path contains two cells in any row or column in which it contains any cells, 
Accordingly, the proposed change of flow rates does not affect the row or column totals except 
for row M and column m, the totals of which it decreases by q. This change therefore amounts 
to using an additional production rate q at the plant A to raise the consumption rate at the mar- 
ket C by q. The amount q is restricted by the condition that the change must not lead to a nega- 
tive flow rate in any cell. Since the flow rates are decreased in the cells Cb and Mc, q is given 
by the lowest flow rate in these "negative cells," i.e., q = 3. 

The program resulting from this change is shown in Table 2; it represents the solution 
to the modified problem for which the consumption rates at a and c are 9 and 8, respectively. 
It is readily verified that no further increase in these consumption rates can be achieved by 


another application of Step 1. The time has therefore come to use Step 2. 
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4, STEP 2: RAISING PRICES 

In this step, the prices at selected plants and markets are raised by the same amount in 
such a manner that (1) the flow rates recorded at the end of the preceding step are compatible 
with the new prices, and (2) new routes become available because their price differentials have 
risen to equal their specific shipping costs. The “vailability of new routes may open up the 
possibility for a successful application of Step 1. If this should not be the case, Step 2 will have 
to be applied again, until enough routes have been brought into play to make Step 1 possible. 

With reference to Table 2, the increase in the consumption rates of a and c that will 
ultimately be achieved by this Step 1 must be furnished by plant A or the store S. To encourage 
increased use of the productive capacity of A, the price at this plant will be kept fixed. Of 
course, the prices at the store and dump are fixed by the nature of the problem. Check marks 
are placed at the ends of rows A and S and at the bottom of column m to indicate the constancy 
of these prices. 

Considering rows A and S against which check marks have been placed, we now use 
condition (1) above to find columns whose prices cannot be allowed to rise. Aside from the 
cell in the dump column, the only cells of row A that contain circles are Ab and Ad. If the 
price of column b were allowed to rise, the price differential for all Ab would exceed the 
specific shipping cost, in spite of the fact that this cell does not carry capacity flow. Condition 
(1) requires therefore that the price of column b be kept constant. We record this by a check 
mark placed below column b. Since it will be found convenient to exhibit the order in which 
check marks are placed against rows and columns, we place this check mark at a lower level 
than the mark already entered below column m. Note that, as far as cell Ad is concerned, the 


price of column d could be raised because this cell carries capacity flow. Applying the same 
argument to row S, however, we find that the price of column d must be kept constant. This is 
recorded by a check mark below column d, which is placed at the same level as that below 
column b. 


We now investigate the columns b and d against which check marks have been placed 
during the preceding investigation of rows A and S. Aside from cells in these rows, the only 
cell in column b containing a circle is cell Cb. Since, however, this cell does not contain any 
flow, a raise of the price of row C would not violate condition (1). Investigating column d in 
this manner, we find that the price of row B must be kept constant, because the price differen- 
tial for cell Bd would otherwise sink below the specific shipping cost in spite of the fact that 
this cell carries a positive flow. We record the constancy of the price of row B by a check 
mark at the end of this row, which is placed to the right of the marks previously made against 
rows A and S. 

We finally investigate row B against which a check mark has been placed during the 
preceding investigation of columns. Aside from cells in the already checked columns b, d, and 
m, the only cells in row B containing circles are Ba and Bc. Since, however, these cells carry 
Capacity flows, condition (1) does not exclude price raises for columns a andc. As no new 
check marks have been made during this second investigation of rows, the checking process 
has reached its end. Our marks indicate that the prices of row C and columns a and c may be 
taised by the same amount p without violating condition (1). 

Next, we have to determine the minimum value of p which will make new routes avail- 
able. Since the prices of rows A, B, and S are kept fixed while those of columns a and c are 
raised, we must watch the cells at the six intersections of these rows and columns. Consider, 
for instance, cell Aa. For this cell to come into play, the price of column a would have to be 
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raised by 2. Checking all six cells in this manner, we find that the minimum value of p that 
will introduce a new route is p= 1. This price increase introduces the cell Sa. The new prices 
are recorded in Table 3, and cell Sa is marked by a circle. 
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MARKETS 


b 
(3) 





















































Before investigating the possibility of a successful Step 1, we must check whether the 
price raises eliminate any cell previously marked by a circle but carrying no flow. The only 
such cell on Table 2 is cell Cb, and the new prices are seen to preclude a positive flow along 
this route. Thus, cell Cb in Table 3 is no longer marked by a circle. 

Similarly, we must check whether or not the new prices require a relabelling of any 
cell that previously carried capacity flow but had a price differential equal to the shipping cost, 
i.e., a cell containing a circle with a slanting line. As a consequence of a price change ina 
column, the price differential of such a cell could rise above the specific shipping cost, so that 
the circle in the cell would have to be marked by a cross. No changes of this type are required 
in our example. 

We are now ready to try Step 1. It is seen that the only possible step of this kind starts 
in cell Sm and brings the flow rate 2 to cell Sa. The resulting program is shown in Table 3; it 
represents the solution to the modified program with the consumption rate 8 at the market c. 
As this consumption rate cannot be raised by a further application of Step 1, another Step 2 is 
due. The appropriate checkmarks are shown in Table 3. The minimum price raise that intro- 
duces a new route is found to be equal 2. The program resulting from this price raise and a 
subsequent Step 1 is shown in Table 4, which represents the solution for the considered numeri- 
cal example, because all entries in row M have now been reduced to zero. 


5. STEP 0: GETTING STARTED 


Steps 1 and 2 have been explained in connection with advanced phases of the process. 
Table 5 shows how the process is started. 
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TABLE 4 
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The price of any row is set at zero, and the price of any column is set equal to the 
smallest shipping cost recorded in this column. In this way no cell is given a price differential 
exceeding the specific shipping cost. Circles are now inserted into the cells where the price 
differential equals the shipping cost. Beginning with row A and proceeding from left to right, 

We now distribute the productive capacity of each plant into the circled cells of the row, paying 
ineach case attention to the limits set, on one hand, by the capacity of the cell and, on the other 
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hand, by the consumption rate of the column and the flow rates already inserted into cells of 
the column. Whenever a cell is given capacity flow, its circle is marked by a slanting line, 
Table 5 shows the program obtained in this manner; it represents the solution to the modified 
problem for which the markets a, c, and d have the consumption rates 9, 5, and 4, respectively, 
The consumption rates at these markets must then be raised to the prescribed values 
by repeated application of Steps 1 and 2. The reader will readily verify that the program of 
Table 5 cannot be improved by Step 1, and that the first application of Step 2 raises the prices 
of row C and columns 4, b, c, and d by 1, bringing cell Bd into play. This price raise requires 
that the circles in cells Ad, Ba, and Bc be crossed and the circle in cell Cm be erased. The 
subsequent Step 1 transfers the flow rate 2 from cell Bm to cell Bd, thus reducing the entry in 
cell Md by 2, i.e., increasing the consumption rate at the market d by 2. Since no further im- 
provement can be made by applying Step 1, prices must be raised once more. It is found that 
the prices of rows B and C and columns a, b, c, d must be raised by 1. This brings cells Ab 
and Sd into play and requires the elimination of the circle in cell Bm. The obvious Step 1 now 
consists in giving the cell Sd the flow rate 2, thus reducing the entry in the cell Md and reach- 
ing the prescribed consumption rate at the market d. The resulting program is that of Table 1. 


6. THE CHECKING PROCESS 

Whereas we have described a systematic way of determining the rows and columns 
whose prices must be kept constant during Step 2, we did not give a method of finding the cells 
affected by Step 1. It turns out that the checking process of Section 4 can also be used for this 
purpose. 

To show this, let us assume that we had overlooked the possibility of applying Step 1 to 
Table 1 and had instead tried to apply Step 2. The resulting check marks are shown in Table 1. 
We note that there appears a check mark at the bottom of column c, the consumption rate of 
which is still below the prescribed value. This indicates that application of Step 2 would be 
premature, because there must be a path of cells that is suitable for Step 1 and terminates in 
cell Mc. Indeed, a check mark was placed at the bottom of column c because this column con- 
tained a circled cell (Cc) with or without slanting line in a row (C) that was checked during the 
preceding investigation of rows. This row was checked because it contained a circled cell 
(Cb) with or without slanting line in a column (b) that was checked during the preceding inves- 
tigation of columns. Continuing in this fashion, we find the cells (Am, Ab, Cb, Cc, Mc) that are 
to be used in Step 1. 

Note that no check marks would have been placed against row C and column c if cell Cb 
had contained a crossed circle or an empty circle and no flow. In the first case, the price dif- 
ferential for this cell would have exceeded the specific shipping cost. As Step 1 requires a 
decrease in the flow rate of this cell, the new shipping program would not be compatible with 
the old prices. In the second case, this decrease would lead to a negative flow rate. Thus, 
circumstances that eliminate the check mark against column c also eliminate the possibility of 
a successful Step 1 terminating in this column. 

It is worth noting that the checking process need not be completed but may be stopped 
as soon as it produces a check mark against a column for which the prescribed consumption 
rate has not yet been reached. A path of cells that terminates in the M cell of this column and 
is suitable for Step 1 can then be found in the manner just described. Thus, the checking 
process either indicates that Step 1 can be applied or furnishes the necessary information for 
the price raises in Step 2. As each of these steps lets the consumption rates at the markets 
approach the prescribed values, a solution must eventually be obtained. 
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1, CONCLUDING REMARKS 

The reader familiar with the method of Ford and Fulkerson [3] or Boldyreff [4] will 
have recognized that the method described here is essentially identical with theirs. The 
weights used by these authors can be identified with the negative prices at the production cen- 
ters and the prices at the consumption centers. Our price differential for a route therefore 
corresponds to their sum of the weights at the endpoints of this route. While this is a simple 
mathematical transformation, it obscures the practical significance of the various steps. Aside 
from bringing out this significance, the present approach has the advantage of applying to the 
generalized transportation problem in which the prices are prescribed at some centers, while 
the production or consumption rates are specified at the remaining centers. 
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